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STATUS  OF  PROGRAM  | 

PHASE 

Origin 

•  TASK 

Operational 

C.  HIT  MOO*  ;  - 

The  program  is  written  in  G635  time-share  series,  FORTRAN  IV,  and  is 
part  of  a  Conversationally  Oriented  Real-Time  Program-Generating  System 
(CORPS).  The  program  consists  of  a  main  program,  three  function  subprograms, 
and  six  subroutines.  The  specific  function  and  relation  of  all  these  routines 
can  be  found  in  the  FUNCTIONAL  FLOW  CHART,  PART  II:  COMPUTER  FUNCTIONAL 
DESCRIPTION  of  this  abstract. 


D.  EQUIPMENT  DETAIL* 

The  program  was  developed  and  is  operational  on  the  WES  G635,  Vicksburg,  MS . 
It  is  also  operational  on  HIS66/80,  Macon,  GA ,  and  Boeing  CDC,  Seattle,  WA. 


i:  iMPkiT-aofNif - 

All  inputs  for  program  H5320  are  either  read  from  a  permanent  data  file  or  are 
cued  and  read  upon  entry  to  the  input  subroutine  DATA  F.  All  specific  input/ 
output  requirements  for  program  H5320  are  given  in  PART  II:  COMPUTER 
FUNCTIONAL  DESCRIPTION  of  this  abstract. 


P.  ADDITIONAL  REMANKI 

Complete  documentation  of  this  program  is  available  from  the  Engineer  Computer 
Programs  Library,  Technical  Information  Center,  WES. 
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B.  PROGRAM  SPECIFICATIONS: 


inguage:  ANSI  FORTRAN  (FORTRAN  IV) 


Solution  Requirements:  The  run  command 


RUN  WESLIB/CORPS/H5320.R 


plus  the  necessary  input  data. 

Method  of  Analysis:  A  finite  difference  scheme  is  used  to  determine 
the  flow  rate  through  the  culverts  and  water  surface  elevation  in  the 
chamber  during  a  lock  operation. 

Core  Requirements  G635 :  21  K  words. 

External  Storage:  Type  of  storage  -  disk.  A  user  may  have  an  input 
permanent  data  file  with  a  maximum  of  1  K  G635  ASCII  words.  If  the 
user  does  not  have  an  input  data  file,  one  of  the  same  size  and  user 
named  will  be  created  during  a  run  of  this  program.  If  a  user  de¬ 
cides  to  save  output  for  graphics  and/or  other  use,  a  user  named 
output  file  with  a  maximum  of  13  K  G635  ASCII  words  will  be  created. 
If  the  output  file  exists,  then  the  existing  output  file  is  used. 
General  Equation:  The  governing  equation  is: 

(k!  +  k2  +  kv(t)  +  k3  +  V  =  etZref  "  z(t)]  "  g  (dt  ) 


Where 


k^  ,  kg  ,  k^  ,  and  k^  =  constant  loss  coefficient  values 


for  the  intake,  upstream  conduit,  downstreaA copduit ,  and. exit 
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k^Ct)  =  the  valve  loss  coefficient  which  is  a  variable  during 
valve  opening  and  has  a  preset  value  of  0.1  for  a  full 
open  valve 

V(t)  =  velocity  of  flow  at  a  reference  location  in  the 
culvert(s)  at  time  t  ,  fps 

p 

g  =  acceleration  of  gravity,  32.2  ft/sec 

e  =  operation  type  (e  =  1  for  filling;  e  =  -1  for  emptying) 
Zref.  =  upper  pool  elevation  (filling)  and  lower  pool 
elevation  Z  (emptying),  ft-datum 

Jj 

Z(t)  =  water-surface  elevation  in  lock  chamber  at  time  t  , 
ft-datum 

L  =  culvert  length,  ft 
dv 

—  =  rate  of  change  of  velocity  of  water  in  the  culvert  at 
at 

2 

time  t  ,  ft /sec 

Range  of  Quantities:  Unlimited  for  practical  applications. 
Accuracy:  Governed  by  accuracy  of  input  data;  the  accumulative 
velocity  error  is  a  measure  of  the  accuracy  of  the  solution. 
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REF :  ER  1110-1-10  -  ENGINEERING  AND  DESIGN  -  Engineering  Computer  Program 
Library  Standards  and  Documentation,  Appendix  B 

PART  I:  ENGINEERING  DESCRIPTION 

1.  PROGRAM  NUMBER:  722-F3-RG-5J1 

2.  TITLE:  H53 20  -  Lock  Filling  and  Emptying  Symmetrical  Systems 

3.  REVISION  LOG:  N/A 

•  PURPOSE  OF  PROGRAM:  The  purpose  of  the  program  is  to  determine  the 
hydraulic  performance  of  a  lock  culvert  system.  The  core  of  the  program 
is  a  finite  difference  scheme  which  is  used  to  determine  flow  rate 
through  the  culverts  and  water-surface  elevation  in  the  chamber  during  a 
lock  operation.  Once  these  are  known,  other  important  hydraulic  param¬ 
eters  are  evaluated.  The  program  output  contains  the  following  computed 
quantities. 

a.  Single  Values. 

(1)  Operation  time  (T),  min. 

(2)  Overtravel  (d  ),  ft,  and  time  (T  ),  min.,  of  occurrence. 

o  m 

b .  Time  Dependent  Values. 

(3)  Time  (t),  min. 

(h)  Valve  opening  ratio  ( b | B ) . 

(5)  Total  inflow  (Q),  cfs. 

(6)  Velocity  head  (h  ),  ft,  at  the  reference  area. 

(7)  Velocity  head  (h  ),  ft,  just  downstream  from  the  valve. 

(8)  Water-surface  elevation  (z),  ft-datum,  in  the  lock  chamber. 

(9)  Water-surface  elevation  (z  ),  ft-datum,  in  the  upstream 
valve  well. 
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(10)  Water-surface  elevation  (z^),  ft-datum,  in  the  down¬ 
stream  valve  well. 

(11)  Piezometric  head  (h  ),  ft-datum,  at  the  culvert  roof  just 
downstream  from  thePvalve. 

(12)  Cavitation  index  (a). 

(13)  Average  velocity  (V),  fps,  in  the  culvert  at  the  reference 
area. 

(lit)  Rate-of-rise  (d  /d^),  ft/min. ,  of  the  lock  chamber  water 
surface. 

(15)  Head  losses  (H£q,  H H£V,  Hj^j,  H ji  1+ »  and  H^qi),  ft,  for 
culvert  intake,  upstream  conduit,  valve,  downstream  con¬ 
duit,  culvert  exit,  and  overall  culvert  system(s), 
respectively. 

(16)  Effective  inertia  (H  ),  ft,  in  the  flow  in  the  culvert. 

0 

(IT)  Overall  drop  (Ht  =  Hjq  +  +  %v  +  +  He),  ft, 

in  water-surface  elevation. 

5.  STEP  SOLUTION: 

a.  Mathematical  Formulation.  The  distribution  of  total  head  at  an 
instant  during  a  lock  filling  operation  is  assumed  to  be  as  shown  in 
Figure  1;  a  comparable  distribution  is  assumed  for  an  emptying  opera¬ 
tion.  The  manner  in  which  each  item  shown  in  Figure  1  is  used  in  set¬ 
ting  up  the  mathematical  model  is  described  in  Table  1.  Two  simulta¬ 
neous  first-order  differential  equations  are  solved  as  a  function  of 
time  (t);  these  equations  are  either  equations  (la)  and  (2),  or  (lb) 
and  (2),  in  the  following  listing. 

(l)  Conservation  of  Energy. 

KV2 

(a)  Normal  Flow  —  =  e(Zrgf  -  z)  -  He  (la) 
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or 

-fcv2 

(b)  Reverse  Flow  — —  =  e(Z  „  -  z)  -  H  (lb) 

2g  ref  e 

(2)  Continuity  (Normal  and  Reverse  Flows) 

f  - e  V  ^ 

The  variables  in  equations  (la),  (lb),  and  (2)  are  as  listed  in  Table  1 
with  the  following  exceptions: 

K  =  k,  +  k  +  k  +  k_  +  k,  ;  where  and  K  are  known  functions 

1  2  v  3  4  v 

of  time. 

e  =  1  for  filling 

-1  for  emptying 

A  =  A/n  A  ;  where  AT  is  the  surface  area  in  the  lock  chamber 
r  L  c  c  L 

and 

dz/dt  =  the  rate-of-change  of  the  elevation  of  the  chamber  water 
surface;  dz/dt  is  positive  for  a  rising  water  level. 

The  only  boundary  conditions  prescribed  are  those  when  t  =  0, 

z  =  Z  for  filling 

Jj 

for  emptying 

V  =  0 
K  =  °° 

b.  Time  Steps.  A  known  time  (t^)  is  advanced  by  a  constant  value 
(At)  to  a  new  time  =  t^  +  At).  Values  of  i  that  are  of  specific 

interest  are  listed  below. 

Value  of  i  Condition 

1  Initial  value  (i.e.,  t^  =  0.0) 
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Value  of  i 
N 

Ml 

M2 

M 

Known :  V .  , 

-  a. 

1st  Prediction 

2nd  Prediction: 

Corrector : 

Letting  superscript 


Condition 

Valve  is  fully  open  (i.e.,  k  =  for  i  _>  N) 

Lock  operation  is  complete 

Filling:  ^  <  Zu  <  zm 

Emptying:  >  ZL  > 

Maximum  overfill  or  over empty 

Filling:  ZM2_1  <  >  ZM2+1 

Emptying:  >  Z^  < 

Final  value  (i.e.,  t  =  (M  -  1)  At) 

M 

z.  ,  K.  .  ,  K.  _ 

1  l+l  *  1+2 


Vi  *  W  V  Vi’ 


zi+1  from  equation  (h) 


Vi+2  Fp^Zi+l’  Vi+1’  Ki+2^ 


At 


Iteration  loop  A 


V.A1  =  F  (z.,  V.  and  V_„,  K.  ,,  ) 
l+l  p  i  i  i+2  l+l 


z^+1  from  equation  (li) 


j  =  1,2,3,  .  represents  the  number  trials 


(loop  A  above),  the  sequence  is  terminated  when 
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(Vi)J  -  (zi+i)J-1  ^1(1°4) 

In  this  program  M  is  the  overall  control.  For  example,  if  M  is  pre¬ 
set  to  a  value  less  than  N  the  solution  will  terminate  before  the  valve 
is  full  open  (i.e.,  M  <  N);  filling  time  and  overtravel  will  not  be 
evaluated. 

c.  Numerical  Integration  (Predictor-Corrector) .  Equations  (2a)  and 
(2b)  are  considered  in  the  form 


V2 

i+1 


Ki+i  L 


In  \  L  fly 

e  Zref  zi+l  g  dt . 


i+1 


Equation  (2),  using  a  first-order  difference,  is  approximated  by 


V  +  V. 

z.  ,  =  z.  +  e  — — - -  A  At 

l+l  l  2  r 


(3) 


(U) 


The  predictor  is  based  on  the  backwards  difference 


dV 


dti+l 


V.  -  V. 

l+l  l 
At 


(5) 


Whereas  the  corrector  uses  the  centered  difference 


dV 


dti+l 


V  -  V 
i+2  i 

2At 


(6) 


Substituting  (k)  and  (5)  in  (3),  and  (U)  and  (6)  in  (3),  yields  the 
predictor,  F^(z,V,K)  and  the  corrector  functions  Fc(z,V,K)  , 
respectively.  The  algebraic  relationships  are  shown  in  Table  1;  the 
sequence  of  calculations  at  each  time  step  is  as  follows. 
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KA V I V I 

10.  Exit  loss,  ■  — - -  Filling  Loss  at  the  ports  into  the  lock 

8  chamber,  is  constant 

Emptying  Loss  at  the  downstream  exit  from 
the  culvert;  is  constant 

11.  Inertia  "head"  Both  The  inertia  "head"  is  placed  at  the 

8  downstream  end  of  the  culvert  (i.e., 

along  with  H^);  L  is  an  effective 
culvert  length,  g  is  the 

acceleration  of  gravity,  and  dv/dt 
is  the  rate-of-change  of  velocity, 

V,  with  respect  to  time,  t. 


i 


I 

*  , 


; 


i 


53; 


gw 


NSTART=0 


nstart=i 


ENTRY  MAIN1 


00  TO  depending  IDEC 


PRINT:  no  perm  file 
re-enter  name 


PRINT:  file  prob 
call  FTS  542-2U03 


READ  INPUTS  from  perm  file:  type  of  operation,  no.  culverts,  title,  upper 
pool  WS  elevation,  lower  pool  WS  elevation,  culvert  roof  elevation  at 
valve,  culvert  length,  width  and  height,  chamber  length  and  width,  loss 
coefficient  for  intake,  US  conduit,  DS  conduit,  and  exit,  time  step,  no 
calc  steps  desired,  no.  of  valve  opening  ratios,  valving  option  and  output 
parameters 


j READ  INPUTS  from  I 

1  perm  data  file :  l  T  / 

valving  \ 

•V,  F. 

I  valve  opening  r  \ 

point  by 

/  ' 

ratios  / 

’^sDoint^' 

c - 

READ  INPUTS  from  perm  data  file:  valve 
sag  time,  first  valve  rate,  time  valve 
stops  motion,  duration  of  stop,  final 
valve  ratio 


NSTART 


Change  \ 
any  data  before 
start  _ 


RUN1,  secondary 
ENTRY  to  DATAE 


return 
^  =  1  ^ 


F 


F 


S 


set  initial  conditions  shift  6 
bits  of  sense  sw  word  ISTAT, 
starting  with  bit  6,  to  far  right 
and  zero  all  other  positions; 

IDEC  =  3 

velocity  prediction  function  use 
backwards  difference  to  compute 
velocity  in  culvert  at  a  given 
time 

_ 

SV  ^ 

-  octal  it 
^or  octal  37 


file  exist  or  is 
already  in  AFT, 
set  IDEC  =  1 


file  does 
not  exist, 
IDEC  =  2 


/SV  =  \ 

octal  05 


compute  water  surf  elev  in  back 
chamber  using  velocity  at  a 
given  time 


file  busy 
IDEC  =  U 


octal  0U 


illegal 
char  in 
file  name 
IDEC  =  5 


/  sv  =  ' 

social  3>* 


velocity  correction  function  use 
central  difference  to  compute 
culvert  velocity  at  a  given  time 


WRITE  to  perm  data  file:  a  header  to  contain  the 
format  of  hov  time  related  variables  are  written, 
number  of  variables  and  what  each  variable  is. 

This  info  will  be  read  by  CORPS  graphic  routine  or 
by  use  with  his  own  routine 


Compute  valve  opening  ratios  and 
number  of  time  steps  required  to 
fully  open  valve 
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3.  EQUIPMENT  AMD  OPERATING  SYSTEM:  The  program  was  developed  on  a  G635 
time-share  system  in  which  input/output  equipment  consisted  of  a  Model  33 
remote  teletype.  It  is  now  operational  on  the  WES  G635>  Vicksburg,  MS; 
HIS  66/80,  Macon,  GA;  and  Boeing  CDC,  Seattle,  WA. 

U.  INPUT  REQUIREMENTS:  The  required  inputs  are  entered  via  the  user's 
time-share  terminal  device  in  free  field  format  or  are  read  from  a  per¬ 
manent  input  data  file.  If  the  inputs  are  entered  through  a  terminal 
device,  then  the  input  subroutine  DATAF  will  supply  all  necessary  cues 
and  reads  for  the  input  data.  A  user  named  permanent  file  will  be 
created  by  DATAF  and  the  inputs  written  to  it.  All  the  inputs  from 
DATAF  are  passed  to  the  computation  subroutine  H5320  via  the  COMMON 
statement.  If  the  inputs  are  entered  from  an  existing  data  file,  then 
the  inputs  are  read  directly  by  H5320.  H5320  also  handles  all  the 

required  output. 

5.  SECONDARY  STORAGE  INPUT  FORMAT: 

a.  The  following  formats  are  used  for  the  input  data  file  whether 
it  already  existed  or  was  created  by  a  run  of  program  H5320 . 

To  write: 

56  FORMAT(I3,2F4.0,A35/I3,3F8.3/I3,5F8.3/I3,Uf8.3/I3,F5.1.2IU,5I2) 

57  FORMAT (I 3, 10 F6, 3) 

To  read: 

10100  FORMAT ( 3X,2F4. 0 ,A35 ) 

10090  FORMAT(3X,3F8.3/3X,5F8.3/3X,UF8.3/3X,F5.1,2lU,5I2) 

10110  FORMAT(3X,10F6.3) 
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Refer  to  line  numbers  2Ul60-2bl80  and  21+200  or  21+220  for  variable 
matchup  of  the  above  write  formats.  Refer  to  20690-20710  and  20710  or 
20760  for  read  formats.  These  line  numbers  are  in  the  source  listings 
on  pages  50,  1+3,  and  1+1+ .  An  example  data  file  is: 


100  1.  2.  EXAMPLE  2 

102  262.800  169.500  128.000 

104  390.000  12.000  14.000  635.400  86.000 

106  0.200  0.060  0.  180  0.350 

108  10.0  31  17  2  l  0  1  3 

110  0.120  1.000  0.500  2.000  3.000 


b.  If  the  user  decides  to  save  the  time-related  output  for  graphics 
or  other  use,  then  an  existing  user  data  file  or  a  newly  created  user 
named  file  is  written  to  in  the  following  manner. 

(1)  Subroutine  HEADER  writes  all  of  its  header  information  to 
the  file.  This  header  information  is  used  by  the  CORPS  system  for 
graphics  use.  If  the  user  wishes  to  read  only  the  time-related  output 
from  the  data  file,  use  the  format  as  described  in  the  second  line  of 
header  information.  Refer  to  line  numbers  21+900-25020  of  the  source 
listing  on  page  52  for  a  complete  description  of  header  formats. 

(2)  The  time-related  variables  are  written  by  the  following 

format . 

10U30  FORMAT(I3,3F6.2,f8.3,F9.0,Uf8.3/F9.2,2F6.2,F8.2,5F8.2) 

Refer  to  line  numbers  22100-22120  in  the  source  listing  on  page  1+6  for 
variable  matchup  to  the  format.  An  example  file  is: 
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03  H5320  13  IS  EDGE 

09  <  1  5</),  (  3x,  3F6.  2,  F3 . 3,  F9 . 0,  4f3 . 3/F9. 2,  2f6.  2,  f3.2,  5f8. 2)  > 

10  31  31  31  31  31  31  31  31  31  31  31  '  31  31  31  31 

11  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1 

12  CVP"E  DESIGNATION?  FO°  H5320  APE* 


13  l»TIM£ 

14  2*VAL"E  OPENING  PATIO 

15  3*CONT'?ACTION  COEF. 

16  4«L0G<’»AL"E  LOSS  COEF.) 

17  5=T0tAL  INFLOV 

13  6*L0CK  VATEP  «"PFACE<  V-<?) 

19  7*"/?  ”ALVE  VELL  V-« 

20  3=PIET  HE  EELOV  VALVE 

21  9  =  C/S  ,fAL”E  VELL  V- e 

22  *-NIT<!  FOP  AEO  VE  ’'API  AELEC  APE* 


23 

MIN”T£?= 1 

FT/MI N* 1 1 

24 

NONE* 2,  3,  4,  10 

FT/SEC* 1 2 

25 

CF  ?*  5 

FT* 13,  14,  15 

26 

FT-DA,r,,M=6,  7,8,9 

10* CAVITATION  PAnAMETEP 
I 1-PATE-OF-PI  PE 
1  2=  t»ELOCI  TY  AT  PEF  APEA 
13*"EL0CITY  HT  AT  PEF  APEA 
14«  VELOCITY  HD  AT  ”ENA  CONT 
1  5=HEADL0Sc  TH®»’  valve 
1  6=  TOTAL  HEACLO'P 
1 7-INEPTI AL  TEPM 
18  =  y-e  Cl  FFEPENTI AL 


6,  1  7,  18 


1  31 


i 
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ive 

0. 

0. 

0.30 

4.00 

e. 

169.500 

262.800 

169.500 

1 69. 500 

9.99 

0. 

0. 

0. 

e.  • 

93.30 

93.  30 

0. 

93. 3e 

102 

0.17 

0.11 

0. 73 

2.11 

2190. 

1  69 . 68  6 

261.969 

1 66. 663 

1  70.246 

0.83 

2.  23 

6.  52 

0.66 

95.97 

34.  19 

8  5;  04 

8 . 03 

93.11 

104 

0.33 

0.23 

e.  70 

1.47 

4453. 

1 70. 249 

259. 363 

1 57. 576 

1  72. 568 

0.7  1 

4.  53 

13.25 

2.  73 

104. 51 

79 . 74 

63.  26 

9. 66 

92.55 

106 

0.50 

0.  38 

0.65 

0.93 

7548  . 

1 71.267 

£52.926 

1 33.883 

1  7  7.9  29 

0.39 

7.  63 

22.47 

7.34 

126.83 

75.  53 

8  5.  64 

5.92 

9  1.53 

123 

0.67 

0.  33 

0. 65 

0.93 

7737. 

172.564 

£52.427 

127.374 

179.561 

0.32 

7.33 

23.03 

3.23 

1  33.  29 

79.  35 

69.9  7 

Q.Z'’ 

9  0.24 

112 

0.33 

0.  38 

0.65 

0.93 

7699. 

173.373 

252 . 523 

128. 668 

160.803 

0.  33 

7.6  4 

2  2.9  1 

3.  15 

131.99 

78 . 5-3 

69  .  £9 

-C.  1  7 

86 .9  3 

1 1  2 

1 . 00 

0.  36 

0.65 

0.93 

7644. 

175. 175 

252. 674 

1  30. 596 

13  2.  606 

0.  35 

7.  73 

22.75 

3 .  e4 

130. t l 

77. 46 

87.33 

-e.  20 

6  7.63 

1  1  4 

1.17 

0.  33 

£.65 

0.93 

7533. 

176. 467 

£52.823 

132. 537 

183.  197 

0.37 

7.  72 

22  .  S3 

7.4? 

123.20 

76.  32 

3  6.54 

-0.  2f 

8  6.  33 

1  1  6 

!  .  33 

£.  33 

0.65 

0.93 

7531  . 

177.749 

252. 970 

1  84. 467 

184, 330 

0,40 

7.67 

22.41 

7.30 

1 26. 50 

75.  19 

85.  25 

-0.20 

35.05 

113 

1 . 50 

0. 33 

0.  65 

0.93 

7475. 

1 79.622 

253. 1 1 7 

136.  333 

135.554 

0.4  2 

7. 61 

2  2.25 

7. 63 

124,42 

7  4.07 

o  3. 93 

-0.  £0 

6  3.  78 

122 

1 . 67 

0.  33 

2.  65 

0.93 

7413. 

1  o  fc  .  23  5 

253. 263 

1  3.5 . 28  5 

186.719 

0. 44 

7.55 

22.  C3 

7.  57 

1  22.  55 

72.9  5 

8  2,72 

-  £ .  £  0 

32.51 

122 

1.35 

C .  3b 

0. 65 

0.93 

7302. 

1 6  1 . 5  5:9 

£53. 407 

140. 1 72 

1 o 7,675 

0.46 

7. 49 

21.91 

7.  45 

1  20 . 0:9 

7  1.85 

6  1.40 

-e.  £0 

6  1.20 

124 

2.00 

0.23 

0.65 

0.  93 

7305. 

15  2. 73  3 

253.65 1 

142.046 

189  .£'23 

0.4': 

7 . 44 

2  1.74 

7.  34 

1  18.35 

76 . 75 

3  0.22 

-0.20- 

6  2  .  £  £ 

12  6 

2.17 

6  •  38 

0. 65 

0.98 

7249  . 

154.013 

253. 043 

143,984 

192.101 

2.51 

7.33 

21.57 

7.23 

117.52 

09  .  C6 

78.99 

-0.20 

73 . 78 

123 

0.  33 

0. 33 

0.63 

0.93 

7193. 

136.248 

253.8  24 

1  4  5.  743 

1  9  1  . £y 1 

3. 3  2 

7.32 

21.41 

7.12 

1  1  5. 20 

Co .  53 

77.  7 ( 

-£.33 

7  7.56 

IOC 

2.5  0 

*-  •  Jci 

2  .  6  5 

0.93 

7840. 

13  6.450 

254. 21 2 

1 50. 674 

19 £.244 

f  .  <  " 

7.17 

2  0.95 

6.3? 

1 »e. 35 

65.  04 

7  4.  49 

1.8  2 

70.  35 

122 

2 .  r  7 

U * 

? .  6  6 

0.30 

8  2:0. 

157, 745 

231. 147 

1  5e. 0  1  1 

190.684 

0.39 

^  .  35 

2  4.  48 

9. 25 

1  1  8.  33 

53  .  48 

7  « .  3  3 

4.  67 

7  5.00 

134 

2: .  j  3 

C.  5? 

8.67 

0.6  1 

9645. 

I 39. 256 

246. 68/ 

14  6.  iy/, 

£62. 1 8 1 

0.  53 

9 . 3P 

23 . 72' 

I  2.  7y 

112.  23 

3  1  .  56 

05.00 

5.  48 

7  3,54 

1  36 

3 . 0  0 

e.  56 

2.69 

0.40 

11241. 

1  •  1 . 028 

240.403 

1  43. 7.62 

2 1  0 . 6  J  0 

0..  50 

J  1  . 

33. 4  3 

1  7. 38 

1  1  4. 53 

4  8.57 

05.99 

5.0? 

7  1.77 

123 

3.17 

r .  6  o 

0.  72 

0.  18 

1  5  8  7  3, 

1  4  3.074 

234.  0o  1 

144.x 31 

212.447 

0.31 

13.  If' 

38. 31 

2  2.  79 

I 12. £4 

34.  74 

64.14 

5 . 6  7 

6  4.73 

140 

3 . 3  8 

2.  7 £ 

0.7  4 

-0.04 

1  '»  J-,  5  » 

l  y  5 . 3  1  < 

220.7  c': 

1  32 , 7  85 

21 7 . 53  7 

6.  60 

1A.65 

4  8.3 

£8.47 

104.66 

25.5  4 

C  2 .  5  7 

4.  4  7 

0  7 . 4  1 

142 

3. 50 

J; .  7  7 

0.  78 

-0.27 

l 5629  . 

1  >  7.9  32: 

220. 467 

161.310 

£  0  •  4/5 

£• .  7  7 

15.91 

4  o  .  58 

55.(0 

9  5.75 

17  .  c  5 

61.  19 

:> . « i 

6  4.87 

144 

3 . 6  7 

e .  5  5 

0  .  ri :  . 

-4.51 

1  050  1  . 

2-- 6 . 655 

215.  01  3 

1 74. B  77 

£  2  2. 4v 2 

l.Ci 

l  6  .  J  P 

49  .  1  1 

37.  45 

78  .  !  9 

1  1 . 47 

59 . 6  f 

2.  4  6 

6  2.14 

146 

3.33 

0 . 8  6 

-0.  76 

1  69  4  ?. 

2  2  3 . 49  9 

212.  7  05 

1  7  '..141 

£37. c l 4 

'  c  3 

1  7 . 3(5 

30. 57 

34.7  1 

63.  33 

6.97 

53.20 

1.26 

54  .  32 

164 

4  •  i* 

1  .  fcfc 

0 . 9  ('. 

-1.00 

1  72  e  I . 

£06. 399 

£11.  o: ;: 

£21.976 

2  40.4',-: 

2 . ;  o 

17.51 

51.14 

46.  70 

50.  2  4 

4.07 

5  0'.  57 

-t.13 

5  o .  2 

150 

4.17 

1  .  VS 

1/  •  y  2' 

-i.ee 

1  ('42  2. 

2.::  4 .  2;  ?  'i 

2  15.17 .0 

2  ■’  3 .957 

2  .  7  71 

2. 27 

1  V  .  £2 

50.  36 

39. 33 

43. 62 

3.44 

54.75 

■  1.23 

0  3.51 

1  5? 

4.8  2 

I  .  t’2 

8 .  y  *■ 

-i.ee 

1  655  6. 

2  12.  1  3 1 

21  5.  51", 

£  0 . 044 

£44.  ?  ■  o  7 

2. 44 

\  f>,  IS 'd 

47.17 

3  7 .  54 

46.  34 

3 .  7  6 

52,  14 

-  1 . 33 

5  £.67 

154 

4.  30 

1 . 

0.  >8 

-1.00 

1  60.  j  7. 

£  1  4." 4 0 

2  1  7.9  49 

£  o  y . 6  9  9 

£45.  1  5,  j 

2.  64 

li .  r* 

4  7.5  5 

35.  6C 

43.9  5 

3 . 5  (. 

49 . 4  o 

-  1  ,  37 

47.92 

150 

4.  07 

1 .  c  ? 

0.90 

-i.ee 

1 5647. 

2  1  7 . 536 

22e. 378 

212.474 

2-46.21 

2.86 

15.93 

46.  57 

33.07 

41 . 57 

3. 87 

4  6.36 

-  1 . 3J 

45.  1 

15  > 

4.32 

l  .  00 

fc.  y£ 

-  1  .  £t 

1  521  4. 

226 . £05 

222. 740 

£  i  5.  3 

£  4  7 .  r :  v 

2  .  1  » 

1  5  ,  Aii 

45.25 

3  1 . 74 

39. 25 

3.  13 

44.  1  4 

-1,60 

4 £  .  09 

16b 

3 . 0  < 

1  .  VS 

1  .  9  v 

-1.00 

1 4766. 

£,.£  .  747 

225. 047 

213.21 4 

£43.21 5 

3 .  80 

15.  er> 

4  3.9  o 

2:9.97 

36. 94 

3.  06 

4  1 . 65 

-1,62 

40.0  5 

18 


H5320 


6.  INPUT  DATA  DESCRIPTION:  The  following  names  are  used  for  the  input 
variables  in  program  H5320. 

AAB  -  8  char  name  of  perm  data  file 

DTS  -  time  step,  sec 

IDEL1  -  output  control  to  valve  full  open;  =  1,  print  every  time 
step  value;  =  2,  print  every  2nd  value;  =  n,  print  every 
nth  value 

IDEL2  -  output  control  from  valve  full  open  to  complete  solution; 

=  1,  print  every  time  step  value;  =  2,  print  every  2nd 
value;  =  n,  print  every  nth  value 

IPRIN1  -  output  control  for  engineering  values;  =  1,  print:  conduit 
and  chamber  geometry,  loss  coefs  and  elevs ;  =  0,  no  print 

IPRIN2  -  output  control  for  head  loss  distributions;  =  1,  print; 

=  0,  no  print 

IVALV  -  type  of  valving;  =  1,  point  by  point;  =  2,  programmed 
functions 

M  -  number  of  calculated  steps  desired 

N  -  number  of  values  of  valve  opening  ratios  (point  by  point) 

RUN  -  35  char  problem  name 

SAG  -  sag  at  one-half  valve  time  of  the  nondimensional  valve 

opening  pattern  (programmed  function) 
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TV1  -  first  valve  rate  (programmed  function),  time  taken  for 

valve  from  closed  to  full  open  position  with  no  stop,  min 
to  full  open 

T1  -  time  valve  stops  initial  motion  (programmed  function),  min; 
no  stop,  T1  =  TV1 

TLAG  -  duration  of  stop  (programmed  function),  min;  no  stop, 

TLAG  =  0 

TV2  -  second  valve  rate  (programmed  function),  time  taken  for 

valve  from  closed  to  full  open  position  as  if  no  stop,  min 
to  full  open.  If  actually  no  stop,  TV2  =  TV1,  never  0. 

TYPE  -  type  of  operation;  =  1,  filling;  =  2,  emptying;  =  0  stop 
program 

XBB  -  valve  opening  ratios  (point  By  point) 

XK1  -  intake  loss  coef 

XK2  -  upstream  conduit  loss  coef 
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XK3  -  downstream  conduit  loss  coef 

XKU  -  exit  loss  coef 

XL1  -  culvert  length,  ft 

XL2  -  culvert  width,  ft 

XL3  -  culvert  height,  ft 

XLU  -  chamber  length,  ft 

XL5  -  chamber  width,  ft 

XN  -  number  of  culverts 

ZL  -  lower  pool  elev,  ft-datum 

ZR  -  culvert  roof  elev  at  valve,  ft-datum 

ZU  -  upper  pool  elev,  ft-datum 

7.  OUTPUT  DATA  DESCRIPTION:  All  output  is  generated  within  subroutine 
H5320.  It  includes  the  given  input  data,  plus  the  output  data  as  de¬ 
scribed  in  the  following  variables. 

BB2  -  valve  opening  at  end  of  first  valve  motion  (programmed 
function),  ft;  Refer  to  figure  given  in  (6)  INPUT  DATA 
DESCRIPTION 

Cl  -  cavitation  index 

2  2 

ERROR  -  accumulated  velocity  error,  ft  /sec 

H  -  overall  head,  ft 

HI  -  effective  inertia,  ft 

HLT  -  total  head  loss,  ft 

HLV  -  head  loss  at  valve,  ft 
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I 

Ml 

M2 

NC 

Q 

RR 

T3 

V 

VH 

VHVC 

XBB 

XCC 

XKC 

XKT 

XKV 

XT 

XTM 

XTX 

Z 

ZOW 


-  count  of  the  number  of  time  steps  at  which  accumulated 
velocity  error  is  printed 

-  number  of  steps  to  fill  or  empty  lock  chamber,  if  input  M 
(number  of  calculated  steps  desired)  is  sufficiently  large 
to  complete  solution 

-  number  of  steps  to  peak  overfill  or  under  empty,  if  input  M 
is  sufficiently  large  to  complete  solution 

-  number  of  culverts 

-  total  inflow,  cfs 

-  rise  of  water  level  in  lock  chamber,  ft/min 

-  valve  opening  time,  min 

-  velocity  in  culvert,  fps 

-  velocity  head  in  culvert,  ft 

-  velocity  head  at  vena  contracta,  ft 

-  valve  opening  ratio 

-  valve  contraction  ratio 

-  equivalent  lock  coef  (based  on  loss  coefs) 

-  sum  of  loss  coefs  (including  valve)  at  valve  full  open 

-  valve  loss  coef 

-  time  at  each  calculation,  min 

-  time  at  extreme  elev,  min 

-  time  to  complete  operation,  min 

-  lock  chamber  water  surface  elev,  ft-datum 

-  water  surface  elev  in  downstream  well,  ft-datum 
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ZOT  -  Lock  chamber  overfill  or  empty,  ft 
ZUW  -  water  surface  elev  in  upstream  well,  ft-datum 

ZV  -  piezometric  head  elev  at  vena  contracta,  ft-datum 

8.  PROGRAM  ERROR  MESSAGES: 

a.  If  file  name  INDATA  is  entered  as  the  name  of  your  permanent 
data  file  and  INDATA  does  not  exist,  then 

NO  FILE  INDATA.  ENTER  FILE  NAME  OR  STOP, 
is  printed  from  subroutine  H5320.  This  is  repeated  until  an  existing 
file  is  attached  or  the  reply  STOP  is  entered.  If  STOP  is  entered,  then 
control  is  returned  to  the  main  program.  Since  program  only  checks  file 
to  see  if  it  exists,  "be  sure"  file  is  your  input  data  file.  If  not  your 
data  file,  program  will  attempt  to  read  bad  data  and/or  abort. 

b.  If  a  program  or  system  error  occurs  in  an  attempt  to  attach  your 
data  file,  then 

FILE  PROBLEM.  CALL  FTS  51+2-21+03  AND  GIVE  THIS  NUMBER  ISTAT 
(where  ISTAT  is  a  12  digit  octal  number) 

is  printed  from  subroutine  H5320.  Control  is  then  returned  to  the  main 

program. 

c.  If  program  cannot  create  or  attach  an  output  file  for  your  re¬ 
quest  of  an  output  data  file,  then 

FILE  PROBLEM; ILLEGAL  CHAR;BUSY;OR  SYSTEM  PROBLEM 
OUTPUT  VALUES  NOT  SAVED  AS  REQUESTED 

is  printed  from  subroutine  H5320.  Control  is  returned  to  the  main 
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d.  If  first  valve  rate  is  less  than  one  time  step,  then 

YOUR  1ST  VALVE  RATE  IS  <  ONE  TIME  STEP 

ENTER  VALVE  OPENING  RATIO  AT  END  OF  1ST  TIME  STEP 

is  printed  from  subroutine  VAL2.  The  ratio  is  read  and  program  continues 

to  next  statement. 

e.  If  program  cannot  create  or  attach  your  input  data  file,  then 
FILE  PROBLEMS; ILLEGAL  CHAR; BUSY; OR  SYSTEM  PROB. 

is  printed  from  subroutine  DATAF.  Control  is  returned  to  the  main 
program. 

9.  VARIABLE  DEFINITIONS: 

a .  Main  program 

HFILE  -  5  char,  name  of  program;  passed  to  WESLIB,  count  routine 
HACCT 

LQZ  -  equal  1,  execute  all  input  cues  and  reads;  equal  2, 

call  WESLIB  routine  RERUN  and  enter  only  desired  inputs; 
initialized  to  1  in  main  program  for  first  run, 

2  thereafter 

LQX  -  equal  1,  print  instructions  for  RERUN;  equal  3  no  print; 

initialized  to  1  in  main  program 
ZZZZZ  -  2  char;  equal  RE,  rerun;  equal  ST,  stop 

b.  Subroutine  H5320,  main  computation  routine 
AAA  -  1  char;  equal  to  Y  for  yes  and  N  for  no  for  response  to 

yes  or  no  question 

AAB  -  8  char,  name  of  permanent  data  file 

2k 


H5320 


AAC 

AAD 

AAE 

AR 

B 

BB1 

BB2 

C 

CFILE 

Cl 

CORR 

D 

DTM 

DTS 

E 

ERROR 

F 

H 


10  char,  used  for  attaching  permanent  data  file 
8  char,  name  of  output  file  for  graphics  and/or  other  use 
10  char,  used  for  attaching  output  file  for  graphics 
and/or  other  use 

ratio  of  total  culvert(s)  cross-section  area  to  lock 
chamber  water  surface  area 

constant  =  2.2,  used  to  determine  valve  loss  coef 
valve  opening  ratio  at  end  of  first  valve  motion 
(programmed  function) 

valve  opening  at  end  of  first  valve  motion  (programmed 
function)  ft 

constant  =  3.2,  used  to  determine  valve  loss  coef 
25  char,  used  for  creating  output  file  for  graphics 
and/or  other  use 
cavitation  index 

velocity  correction  FUNCTION  subprogram  (central  differ¬ 
ence),  fps 

constant  =  .65,  used  to  determine  valve  contraction  coef 
time  step,  min 
time  step,  sec 

constant  =  .9,  used  to  determine  valve  contraction  coef 

2  2 

accumulated  velocity  error,  ft  /sec 

constant  =  .9,  used  to  determine  valve  contraction  coef 
overall  head,  ft 
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HI 

HLT 

HLV 

HL1 

HL2 

HL3 

HLh 

I 

IDEC 


IDEIJ.  - 


IDEL2 


IPRIN1  - 


IPRIN2  - 


IVALV  - 


effective  inertia  term,  ft 

total  head  loss,  ft 

head  loss  at  valve,  ft 

intake  head  loss,  ft 

upstream  conduit  head  loss,  ft 

downstream  conduit  head  loss,  ft 

exit  head  loss,  ft 

count  of  the  number  of  time  steps 

switch  equal  to: 

1  -  file  exist  or  is  in  AFT  2  -  file  nonexistent 
3  -  file  problem  U  -  file  busy 

5  -  illegal  char  in  file  name 

output  control  to  valve  full  open;  *  1,  print  every  time 
step  value;  =2,  print  every  2nd  value;  =  n,  print  every 
nth  value 

output  control  from  full  open  to  completed  solution; 

=  1,  print  every  time  step  value;  =  2,  print  every  2nd 
value;  =  n,  print  every  nth  value 

output  control  for  engineering  input  variables;  =  1, 
print;  =  0,  no  print 

output  control  for  loss  distribution;  =  1,  print;  =  0, 
no  print 

type  of  valving;  =  1,  point  by  point;  =  2,  programmed 
function 
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IVAL1  -  switch  variable  set  to  0  to  be  passed  to  subroutine 
DATAF 

IVAL2  -  switch  variable  set  to  0  to  be  passed  to  subroutine 
DATAF 

KK  -  switch  variable  set  to  0  to  be  passed  to  subroutine 
DATAF 


LINE  -  line  numbers  in  output  file  for  graphics  and/or  other 
use 

M  -  total  number  of  calculated  steps  desired 

Ml  -  number  of  steps  to  fill  or  empty  lock  chamber  if  M  is 

sufficiently  large  to  complete  solution 
M2  -  number  of  steps  to  peak  overfill  or  under  empty  if  M 
is  sufficiently  large  to  complete  solution 
M3  -  switch;  omit  corrector  sequence  for  low  velocity  flows 
N  -  number  of  valve  opening  ratios 

NC  -  number  of  culverts 

NSTART  -  switch;  =  1  if  input  is  first  read  from  permanent  data 
file,  else  =  0 

PRED  -  velocity  prediction  FUNCTION  subprogram  (backward 
difference),  fps 
Q  -  total  inflow,  cfs 

RR  -  rise  of  water  in  lock  chamber,  fpm 

RUN  -  35  char  problem  name 
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SAG  -  sag  at  one-half  valve  time  of  the  nondimensional  valve 
opening  pattern  (programmed  function) 

T1  -  time  valve  stops  initial  motion  (programmed  function), 
min 

TLAG  -  duration  of  stop  (programmed  function),  min 
T3  -  valve  opening  time,  min 

TV1  -  first  valve  rate  (programmed  function),  time  taken  for 
valve  from  closed  to  full  open  with  no  stop,  min  to 
full  open 

TV2  -  second  valve  rate  (programmed  function),  time  taken  for 
valve  from  closed  to  full  open  as  if  no  stop,  min  to 
full  open.  If  actually  no  stop,  TV2  =  TV1,  never  0 
TYPE  -  type  of  operation;  =  1,  fill;  =  -1,  empty;  =  0,  stop 
program 

V  -  velocity  in  culvert,  fps 

VH  -  velocity  head  in  culvert,  ft 

VHVC  -  velocity  head  at  vena  contracta,  ft 

XBB  -  valve  opening  ratios 

XCC  -  valve  contraction  ratios 

XK  -  sum  of  loss  coefs  (including  valve) 

XKC  -  equivalent  lock  coef  (based  on  loss  coefs) 

XKT  -  sum  of  loss  coefs  (including  valve)  at  valve  full  open 
XKV  -  valve  loss  coef 
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XK1  -  intake  loss  coef 

XK2  -  upstream  loss  coef 

XK3  -  downstream  loss  coef 

XKU  -  exist  loss  coef 

XL1  -  culvert  length,  ft 

XL2  -  culvert  width,  ft 

XL3  -  culvert  height,  ft 

XLU  -  chamber  length,  ft 

XL5  -  chamber  width,  ft 

XN  -  number  of  culverts 

XT  -  time  at  each  calculation,  min 

XTM  -  time  at  extreme  elev,  min 

XTX  -  time  to  complete  operation,  min 

YY  -  working  storage 

Z  -  lock  chamber  water  surface  elev,  ft-datum 

ZCHEC  -  lock  chamber  water  surface  elev  at  time  t.,,  ,  ft-datum 

l+l 

ZOW  -  water  surface  elev  in  downstream  well,  ft-datum 
ZNEW  -  lock  chamber  water  surface  FUNCTION  subprogram  to 
evaluate  successive  elevs,  ft-datum 
ZL  -  lower  pool  water  surface  elev,  ft-datum 

ZOT  -  lock  chamber  overfill  or  under  empty,  ft 

ZR  -  culvert  roof  elev  at  valve,  ft-datum 

ZREF  -  upper  pool  water  surface  elev,  filling;  lower  pool  water 


surface  elev,  emptying;  ft-datum 
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ZU  -  upper  pool  water  surface  elev,  ft-datum 

ZUW  -  water  surface  elev  in  upstream  well,  ft-datum 

ZV  -  piezometric  head  elev  at  vena  contracta,  ft-datum 

c.  Function  PRED,  velocity  prediction  function;  uses  backwards 
difference  to  evaluate  the  culvert  velocity  at  a  given  time. 

AR  -  ratio  of  total  culvert(s)  cross-section  area  to  lock 
chamber  water  surface  area 
DTS  -  time  step,  sec 

PRED  -  symbolic  name  of  single  valued  function  used  as  the  return 
for  the  output  value  to  the  calling  program,  fps 
TYPE  -  type  of  operation;  =  1,  filling;  =  -1,  emptying 

VST  -  culvert  velocity  at  time  t^  or  t^+^  ,  fps 

XB  -  working  storage 

XC  -  working  storage 

XKF  -  sum  of  loss  coefs  (including  valve)  at  time  t^+^  or 
ti+2 

XL1  -  culvert  length,  ft 

ZREF  -  upper  pool  water  surf  elev,  filling;  lower  pool  water 

surface  elev,  emptying;  ft-datum 

ZST  -  lock  chamber  water  surface  elev  at  time  t.  or  t.,.,  , 

1  l+l 

ft-datum 

d.  Function  COOR,  velocity  correction  function;  uses  central  dif¬ 
ference  to  evaluate  the  culvert,  velocity  at  a  given  time. 
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AR  -  ratio  of  total  culvert(s)  cross-section  area  to  lock 
chamber  water  surface  area 

CORR  -  symbolic  name  of  single  valued  function  used  as  the  re¬ 
turn  for  the  output  value  to  the  calling  program,  fps 
DTS  -  time  step,  sec 

TYPE  -  type  of  operation;  =  1,  filling;  =  -1,  emptying 
VAH  -  culvert  velocity  at  time  or  t^+£  ’  fps 

VST  -  culvert  velocity  at  time  t^  or  t^  ,  fps 

XB  -  working  storage 

XC  -  working  storage 

XKF  -  sum  of  loss  coefs  (including  valve)  at  ^  or  t^ 

XL1  -  culvert  length,  ft 

ZREF  -  upper  pool  water  surface  elev,  filling;  lower  pool  water 

surface  elev,  emptying;  ft-datum 

ZST  -  lock  chamber  water  surface  elev  at  time  t.  ,  or  t.  , 

l-l  i 

ft-datum 

e.  Function  ZNEW,  evaluate  water  surface  elev  in  lock  chamber  using 
velocity  at  a  given  time. 

AR  -  ratio  of  total  culvert(s)  cross-section  area  to  lock 
chamber  water  surface  area 
DTS  -  time  step,  sec 

TYPE  -  type  of  operation;  =  1,  filling;  =  -1,  emptying 
VF  -  culvert  velocity  at  time  t^+1  ,  fps 
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VST  -  culvert  velocity  at  time  t^  ,  fps 

ZNEW  -  symbolic  name  of  single  valued  function  used  as  the 
return  for  the  output  value  to  the  calling  program, 
ft-datum 

ZST  -  lock  chamber  water  surface  elev  at  time  t^  ,  fps 

f.  Subroutine  TEST,  determine  status  codes  returned  by  file  system 
for  an  attach  operation. 

IDEC  -  decimal  value  for  status  codes;  IDEC  = 

1  -  file  exist  or  is  in  AFT 

2  -  file  nonexistent 

3  -  file  problem 
-  file  busy 

5  -  illegal  char  in  file  name 

ISTAT  -  status  code  (12  digit  octal  number)  returned  by  file 
system;  converted  to  decimal  as  represented  by  IDEC 
ITEST  -  2  digit  octal  number  =  00,  0h ,  05,  31*  or  37 ;  used  as 

test  against  status  code  ISTAT  to  determine  IDEC  value 

g.  Subroutine  VAL2,  called  if  valve  option  is  programmed  functions. 
Computes  valve  opening  ratios  and  number  of  time  steps  to  fully  open 
valve . 

A  -  sag  at  one-half  valve  time  of  the  nondimensional  valve 
opening  pattern 

BB1  -  valve  opening  ratio  at  end  of  1st  time  step 
DTM  -  time  step,  min 

DTS  -  time  step,  sec 

N  -  number  of  values  of  valve  opening  ratios 
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PI  -  constant  =  3. 1U1593 
TLAG  -  duration  of  stop,  min 

TVI  -  first  valve  rate,  time  taken  for  valve  from  closed  to 

full  open  with  no  stop,  min  to  full 
TV2  -  second  valve  rate,  time  taken  for  valve  from  closed  to 
full  open  as  if  no  stop,  min  to  full.  If  actually  no 
stop,  TV2  =  TVI,  never  =  0 
T1  -  time  valve  stops  initial  motion,  min 

T2  -  time  at  which  valve  motion  starts  following  the  partial- 
open  stop,  min 
T2N  -  working  storage 

T2S  -  working  storage 

XBB  -  valve  opening  ratios 

XT  -  time  at  each  calculated  ratio,  min 

h.  Subroutine  DATAF,  supports  inputs  via  teletype  and  rerun  option. 
All  variables  are  as  described  in  subroutine  H5320  except : 

CFILE  -  2U  char,  used  for  creating  a  permanent  input  data  file 
LQZ  -  equal  1,  execute  all  input  cues  ard  reads;  equal  2,  call 
WESLIB  routine  RERUN  and  enter  only  desired  inputs 
IVAL1  and 

IVAL2  -  switch  pair  for  setting  value  of  rerun  option  LQZ 

IVAL1  and  IVAL2  sets  LQZ 
10  2 
0  12 

all  other  conditions,  switch  ignored 
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KK  -  switch  for  rerun  option  LQZ;  KK  =  1  sets  LQZ  =  2;  any 
other  value  of  KK,  switch  ignored 
KKK  -  total  number  of  inputs,  passed  to  RERUN 
JKL  -  direct  return  from  RERUN  to  desired  input  read 
LQX  -  equal  1,  print  instructions  from  RERUN;  equal  3,  no  print 

i.  Subroutine  PRINTT,  provides  a  listing  of  input  variables  at  the 
teletype  during  a  run  of  program  H5320. 

AAA  -  1  char;  =  N,  no  print  listing;  =  Y,  print 

j.  Subroutine  HEADER,  called  if  desire  output  stored  for  graphics 
and/or  other  use.  HEADER  writes  a  heading  to  the  output  file  which  con¬ 
tains  the  format  the  file  was  written  in,  the  number  of  time  related 
variables  written  and  what  each  variable  is,  and  the  type  of  axis  to 
graph  the  data  on.  The  time  related  variables  are  then  written  under 
the  heading. 

M  -  number  of  calculated  steps 

10.  EXAMPLE  CASE:  Determine  the  hydraulic  performance  of  a  lock 
culvert  system  during  a  filling  operation. 


Input  Data: 

Type  operation 

1 

Problem  title 

Example  1 

No.  culverts 

2 

Upper  pool  elev,  ft-datum 

262.80 

Lower  pool  elev,  ft-datum 

169.50 

Culvert  roof  elev,  ft-datum 

128.00 

Culvert  length,  ft 

390.00 

Culvert  width,  ft 

12.00 

Culvert  height,  ft 

lU .  00 

Chamber  length,  ft 

685 .  ^0 
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Chamber  width,  ft 

Intake  loss  coef 

U.  S.  loss  coef 

D.  S.  loss  coef 

Exit  loss  coef 

Time  step,  sec 

No.  of  stepped  values 

Type  of  valving 

Valve  sag  at  1/2  open  time 

First  valve  rate,  min 

Time  valve  stops  motion,  min 

Duration  of  stop,  min 

Final  valve  rate,  min 

Output  control,  engr  inputs 

Output  control,  loss  distribution 

First  output  rate 

Final  output  rate 


86. 00 
0.20 
0.06 
0.18 
0.85 
10.00 
T3 
2 

0.12 
.  00 
M0 
0.00 
H.00 
1 


Out  put : 


EXAMPLE  1 


CON  CPI  T  GEOMETW:  LENGTH  *  390.0  FT.  WIDTH«  -12.0  FT. 

HEIGHT*  14.0  ft.  N»  2-CrLV£cT  OPERATION) 

CHAMBEP  GEOMETRY:  LENGTH  *  685.4  FT.  WIDTH  *  86.0  FT. 

LO«p  CO  EFF  I  C I  ENT'.  INTAKE  0.200 
MPCT°EAM  CONCPIT  0.060 
DOVNfT^EAM  COMDPIT  0.130 
MANIFOLD  0.350 

TOTAL! VALVE  OPEN)  1.390  LOCK  COEF. =0.843 

EL  EVA  T I  ON  ?!  v^EP  POOL  262.30 
LO''EP  oqql  169.50 
CVLVEPT  poof  123.00 


VAL"E  OPEN?  AT  4.00-MIN.  ^ATE  FOR  4.00  MINPTEF 

ACC. "EL-  EP"OP  TO  1=  58  i?  0.00595 
ACC. "EL.  En"0"  To  1=  64  ip  0.00593 


LOCK  FILLF  TQ  "opEP  POOL  IN  9.56  MINPTEP 
E<T"EME  ELE"ATtON  AT  TIME  =  10.50  MINV’EP 
MAX  COMr"rED  0"EPTPA"EL  =  1.57  FT. 
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VALVE  PA °AMETEPS  ELEVATIONS  (  FT-CAT”M) 


TIME  OPEN. 

CONT. 

LOSS  ' 

TO  TAL 

LOCK 

VS  VAL'fE 

PI  EZ  .HE 

DS  VALVE 

<M )  RATIO 

CO  £F. 

COEF.  : 

INFLOW 

CHAMBEP 

WELL 

VALVE 

WELL 

CFS 

CAV.PAP 

0. 

0.00 

0.30 

10000.00 

0. 

169 . 50 

262.60 

169.50 

169.50 

9.99 

0.17 

0.03 

0.S0 

2147. a  l 

556. 

169. 55 

262. 75 

164. 34 

169. 53 

0.34 

0.33 

0.05 

0.  79 

531.46 

1115. 

169 . 69 

262. 53 

163. 62 

169.8  3 

0.83 

0.50 

0.03 

0.79 

232.25 

1  63  1 . 

169.93 

262. 31 

162.38 

170. 26 

0.80 

0.67 

0.  1  1 

0.  73 

127.64 

2258. 

170.26 

261.92 

160. 56 

1  70.3  6 

0.77 

0.33 

0.  14 

0.76 

79.36 

23  48. 

170. 69 

261 . 39 

1  58 . 0  7 

171.64 

0.  72 

1.00 

0.  17 

0.  75 

53.25 

3452. 

1 71 . 23 

260. 74 

154.80 

1  72. 62 

0.67 

1.17 

0.20 

0.73 

37.63 

4075. 

171.87 

259.92 

150.37 

173.81 

0.60 

1.33 

0.23 

0.70 

29.23 

458  5. 

172.60 

259.  1  6 

151.25 

1  75.06 

0.61 

1.50 

0.26 

0.  63 

22.63 

5147. 

1  73. 43 

258. 21 

148. 24 

1 76. 52 

0.57 

1.67 

0.  30 

0.  65 

17.23 

5799. 

174. 35 

256.97 

140.57 

173.29 

0.47 

1  .83 

0.34 

0.65 

13.00 

6554. 

175.40 

255. 36 

140.19 

130. 42 

0.46 

2.00 

0.33 

0.  65 

9.64 

7424. 

1  76. 59 

253.25 

l  33  .  1  0 

183.03 

0.44 

2.17 

0.42 

0.66 

7.  04 

3414. 

177.93 

250.53 

l  34.9  5 

186.21 

0.40 

2.33 

0.47 

0.  67 

S.  06 

9520. 

179. 45 

247.09 

1  3  f .  43 

190.05 

0.36 

2.50 

0.51 

0. 63 

3. 59 

10726. 

131. 17 

242.86 

128.59 

194. 62 

0.33 

2.67 

0.56 

0.69 

2.51 

1  1994. 

18  3.  10 

237.87 

127.25 

199.92 

0.31 

2.83 

0.61 

0.71 

1.73 

13269. 

185. 24 

232.29 

1 23. 43 

205.82 

0.  32 

3.00 

0.67 

0.73 

1 .  18 

14477. 

187.59 

226.  48 

132.81 

212.10 

0.37 

3.17 

0.72 

0.75 

0.80 

15542. 

190.  14 

220.94 

140. 53 

213.38 

0.46 

3.33 

0.77 

0.  73 

0.53 

16400. 

192.85 

216.  19 

151.06 

224. 29 

0.60 

3.50 

0.33 

0.8  1 

0.35 

17016. 

195.  63 

212. 62 

1 63. 36 

229 . 54 

0.8  2 

3.67 

0.89 

0.84 

0.  23 

1 7369. 

193. 60 

210.40 

1 76. 21 

233.95 

1.12 

3.33 

0.94 

0.37 

0.  15 

1 7542. 

201.57 

209. 47 

188. 57 

237.54 

1 .52 

4.  00 

1.00 

8.90 

0.10 

17548. 

204. 54 

209.43 

199.50 

240.54 

2.03 

4.50 

1.00 

0.90 

0.  10 

16364. 

213.  19 

216.39 

207. 75 

244. 50 

2.51 

5.00 

1.00 

0.90 

0.  10 

15035. 

221 . 18 

223. 62 

21 6. 33 

247. 61 

3.20 

5.50 

1 . 00 

0.90 

0.10 

13704. 

228 . 49 

230.25 

224. 19 

250. 45 

4.  1  0 

6.00 

1.00 

0.90 

0.10 

12373. 

235.  13 

236. 27 

231.33 

253. 03 

5.  30 

6.50 

1.00 

0.90 

0.10 

1 1042. 

241 . 09 

241 . 67 

237. 74 

255. 34 

6.96 

7.00 

1.00 

0.90 

0.10 

9711. 

246. 37 

246. 46 

243. 4 1 

257.40 

9 . 36 

7.50 

1 . 00 

0.90 

0.  10 

8  38  0. 

250.97 

250. 63 

248 . 36 

259. 18 

9.99 

8.00 

1.00 

0.96 

0. 10 

7049. 

254.90 

254.  19 

252. 59 

260. 71 

9  .99 

8.50 

1.00 

0.90 

0.  10 

5718. 

253. 15 

257.  1  3 

256. 08 

261 .97 

9.99 

9.00 

1.00 

e.90 

0.  10 

433  7. 

260. 72 

259. 47 

253.8  4 

262.97 

9.99 

9.5B 

1 . 00 

0.90 

0.  10 

3056. 

262. 61 

261 .  18 

260.38 

263. 71 

9.99 

10.00 

1 . 00 

0.96 

0.  10 

1  725. 

263.83 

262. 28 

262. 19 

264. 18 

9  .99 

10.50 

1 . 00 

0.90 

0.10 

394. 

264. 37 

262. 77 

262. 77 

264.39 

9  .99 

11.00 

1 . 00 

0.90 

0.  10 

-833. 

264. 24 

262. 73 

262. 63 

264. 1 6 

9 .99 

11 .56 

1 . 00 

0.90 

0.10 

-  1  446. 

263. 64 

262. 59 

262. 30 

263. 39 

9.99 

12.00 

1 . 60 

« 

0.90 

0.  10 

-1393. 

262.89 

262. 60 

262. 34 

262. 67 

9  .99 
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PATE 

VEL  HDS<F> 

DIFFERENTIAL  HDS(F) 

PI  SE 

VEL. 

REF. 

VENA 

VALVE 

TOTAL  INERTIA 

OVERALL 

<  F/M) 

<  F/S) 

AREA 

CONT. 

LOSS 

0. 

0. 

0. 

0. 

0. 

93. 300 

93. 300 

0. 

93.300 

0.17 

0.566 

1 . 65A 

0.0A2 

98. A A3 

9 l.  139 

91,243 

2.ee9 

93.253 

0.33 

1.135 

3.  318 

0.171 

99. 1 3A 

90.861 

9  1 . 08  2 

2.029 

93.  11  1 

0.50 

1.712 

5. 00A 

0. 389 

100.317 

90.311 

90.313 

2.061 

92.874 

0.67 

2.  299 

6.721 

0.  701 

102.056 

39.532 

90. A37 

2.1  02 

92.540 

0.8  3 

2.899 

3  .  A75 

1.115 

1  0  A,  A  A  1 

88.516 

39.955 

2.  151 

92. 107 

1.00 

3.513 

10. 273 

1 . 639 

107.573 

37. 2SA 

89. 368 

2.212 

9  1.572 

1.17 

A.  1  A8 

12. 123 

2.  28  A 

11 1 .8  A  1 

3  5.  9A0 

88.887 

2.  0A3 

9  0 . 9  3  A 

1.33 

A.  667 

13.6A6 

2.891 

1  !  0.  79A 

8  A.  529 

83 . 259 

1.932 

90. 199 

1 .50 

5.  239 

15.319 

3.  6AA 

1 13. 615 

32.  A75 

37. 175 

2.  183 

89.374 

1 . 67 

5.903 

1 7. 259 

A.  625 

121.023 

79.9  29 

85.895 

2.  536 

38. AA5 

1  .83 

6.  672 

19.507 

5.909 

121.069 

76.330 

8  A. A52 

2.929 

37.397 

2.00 

7.557 

22.  095 

7.  53  1 

122.732 

73. 06A 

82.8A3 

3.351 

86.212 

2.17 

8. 56A 

25. 0A0 

9.  736 

125. 322 

63 . 51  A 

3  1 .074 

3.  778 

8  A .  8  63 

2.33 

9. 69  0 

23.333 

12. A65 

123.075 

63. 039 

79 . 1  69 

A.  1  67 

83. 3A7 

2.50 

10.918 

31.921 

15.822 

1 30. 096 

56.  760 

77.171 

A.  4  60 

8 1 . 630 

2.67 

1  2. 209 

35.697 

19.73  7 

130. A02 

49.613 

75. 1 39 

A.  58 A 

79 . 702 

2.83 

13. 507 

39. A91 

2A. 21 6 

128.069 

Al .838 

73. 127 

A.  475 

77. 560 

3.00 

1  A. 736 

A3. 086 

23.327 

122. A9 2 

33.988 

71. 174 

A.  096 

75. 206 

3.  1  7 

1  5.820 

A6. 255 

33. 222 

113. 6  3  A 

26.  Al 6 

69.273 

3.  A65 

72.660 

3.  33 

16. 693 

A3 .809 

36.992 

102. 1 26 

19.654 

67. 37A 

2.  653 

69.950 

3.50 

1 7. 321 

50. 6AA 

39.826 

39.090 

1  A. 032 

65. Aes 

1.73  3 

67.116 

3.67 

17.701 

51 . 753 

A  1 . 59  0 

75. 779 

9.662 

63.313 

0.9  A3 

64.  19  7 

3.83 

17.856 

52.209 

A2. 326 

63. 230 

6.  A57 

61.053 

0.  28  7 

61 . 234 

A.  00 

17.863 

52.227 

A2. 355 

52. 290 

A.  236 

58 . 8  74 

-0. 538 

58.257 

A.  50 

16.657 

A3. 702 

36.830 

A5. A69 

3.  63  3 

51.194 

-  1 . 58  3 

A9 . 6 1  0 

5.00 

1 5. 30A 

AA. 7A8 

3 1 . 09  2 

38 . 38  6 

3.  109 

A3. 213 

-  1 . 599 

Al .619 

5.50 

13.950 

A0. 736 

25.831 

31.390 

2.  58  3 

35.905 

-  1 . 599 

34. 306 

6.  00 

1  2.  59  5 

36.825 

21.057 

25.996 

2.  106 

29. 269 

-1.599 

27.670 

6.50 

1 l. 2A0 

32.863 

16.770 

20. 70 A 

1.677 

23.310 

-1.599 

21.711 

7.ee 

9.385 

28.982 

12.971 

16.013 

1 . 29  7 

13. 029 

-  1. 599 

1 6. 430 

7.50 

3. 530 

2A.9A0 

9. 659 

l  1 .92A 

0.966 

13.426 

-  1 . 599 

11.826 

8 . 00 

7.  175 

20.9  79 

6.83A 

8  .  A37 

0. 68  3 

9.  A99 

-  1 . 599 

7.900 

3.50 

5.320 

17.017 

A.  A9  7 

5.  552 

0.  A50 

6.  250 

-  1 . 599 

A.  651 

9 .  eo 

A.  A65 

13.056 

2.  6A7 

3.  2  68 

0.  265 

3.  679 

-1.599 

2. 08  0 

9.50 

3.  1  10 

9 . 09  A 

1 . 23  A 

1. 58  6 

0.  1  23 

1.  78  5 

-1.599 

0.186 

10. 00 

1 . 756 

5.  133 

0.  A09 

0,  505 

0.  0 A 1 

0.  569 

-1.599 

-1.031 

10. 50 

0.  A01 

1.171 

0.02  1 

0.  026 

0.  002 

0.  030 

-  1 . 599 

-  1 . 570 

11.00 

-0.3  A3 

-  2 .  A3  0 

0.09  6 

0,096 

0.010 

e.  1 33 

-1.175 

-  1 . AAA 

11.50* 

-  1 . A72 

-A. 30A 

0.  23  8 

0.  28  3 

0 .  0  29 

0.  A00 

-e. 299 

-e.836 

12.00 

-  1 . Al 7 

-A. 1 AA 

0. 267 

0.  267 

0.  027 

0.  371 

2.  6 49 

-0.093 

TO  YOU  NEED  TO  STOPE  TOE  OUTPUT  IN  A  DATA  FILE?Y  OP  II 
=N 


ENTEP.  PEPUN  OP  STOP 
=PEPUN 

rEP'TI  OPTION  PEPMITS  YOU  TO  CHANGE  ANY  OP  ALL  INPUT  UAPIABLES. 
AT  >>>=  CUE  TYnE  IN  THE  TUO  LETTEPSC AA,  AE/ ETC.  >  CO PPESFON DI N G 
TO  THE  UAPIAELES  YOU  WISH  TO  CHANGE.  THEN  AT  NEYT  =  CUE.  ENTEP 
THE  N1JM  EP I  CAL  VALUE.  TO  TERMINATE  DATA  ENTRY.  TYPE  A  CARRIAGE 
PETUPN  AT  >>>  =  OUE. 


'*  >b- 
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>>  > 

=  AB 

AE-ENTER  YO VP  PROBLEM  NAME 
=  EX  AMPLE  2 
>>  > 

=A0 

AC- ENTER  TOTAL  NO.  OF  CALC.  STEPS  DESIREE 
=  31 
>>  > 

=AV 

AV- ENTER  FIRST  VALVE  PATE 
=  1 
>>> 

=A'-r 

AV- EN TER  TIME  VALVE  STOPS  MO TI ONC M I N V TES> 

=  .5 

>>  > 

=A'C 

A'?  -  EN  TER  DURATION  OF  STO  PC  M  I NVTES  ) 

=  2 
>>  > 

*AY 

AV- ENTER  FINAL  VALVE  RATE 
=  3 
>  >  > 

=RA 

rA- ENTER  2nd  OUTPUT  CONTROL* l = DO. 8= DO  NOT  t^IMT 
“  & 

>>  > 


E\'A-::rLE  2 


CO  N  r  V 1  T  G  EO  M  E  TRY  s  L  EN  G  TO 
H El  CRT* 


CHAMFER  GEOMETRY:  LENGTH 


=  39  2.0  FT.  WIDTH*  12.0  FT. 

lA.e  ft.  fj=  2-C’XVE,>T  OPERATION) 

=  63  5.  A  FT.  WIDTH  =  86.0  FT. 


LOSS  COEFFICIENTS  IN  TATE 
iipFTREAM  CONDUIT 
DO VN S T  V EAM  CONDUIT 
MANIFOLD 
TO  TALC  VALVE  OPEN) 


0.200 
C.0  60 
0 .  18  0 
E.3  5K 

1.390  LOCH  COEF .  =  0. 8  A3 


IT. E”ATI  OM S*  Ut:r,ER  POOL  262.38 
1  OVER  T'OOL  16V.  50 
CUL  ’ERT  ROOF  123.00 


VALUE  ORF.NS  AT  I  •  80-MI  M.  RATI:  FOR  0.80  MI  NVTES 

VAL’LE  HOLDS  FOP  R.  68-MI H.  AT  5.  32ft.  J 
THEN  OPENS  AT  3.00-tJIN.  RATE  TO  FULL  OPEN 
TOTAL  VALVE  TIME  =  4.17-M1N.  <N=  26) 


ACC.VLL.  ERROR  TO  1=  31  IS  0.01  A 08 

ACC • ”EL .  ERROR  TO  1=  31  IS  8.81A00 

LO  C  <  OPERATION  IS  NOT  COMPLETE 
LOCH  rOEr  NOT  COMPLETE  OVERTRATOX 
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VALVE  PAPAMETEPS  ELEVATIONS  (FT-EATm) 


TIME  OPEN. 

CONT. 

LOSS  TOTAL 

LOCK 

US  VALVE 

PI  Et.HD 

DS  VALVE 

<M)  RATIO 

COEF. 

COEF.  INFLOW 

CHAMBER 

WELL 

VALVE 

WELL 

CFS 

"  * 

CAV.PAP. 

0. 

0. 

0.8  0 

10000.00 

0. 

169.50 

262.30 

169. 50 

1 69 . 50 

9.99 

0.17 

0.  1  1 

0. 73 

127.64 

2190. 

1 69 . 69 

261.97 

166.66 

170.25 

0.88 

0.33 

0.  23 

0.  70 

29.23 

4453. 

170.25 

259. 36 

157. 58 

172.57 

0.71 

0.50 

0.33 

0.65 

9.64 

7548. 

171.27 

252.93 

133.83 

177.93 

0.39 

0.67 

0.38 

0.65 

9.64 

7737. 

172. 56 

252.43 

127. 37 

1 79. S6 

0.32 

0.83 

0.33 

0.65 

9.64 

7699. 

173.87 

252.53 

1  28 . 69 

180.80 

0.33 

1.00 

0.38 

0.65 

9.64 

7644. 

175.17 

252. 67 

130.60 

182.01 

0.  35 

1.17 

0.  38 

0.65 

9.64 

7583. 

176.47 

252.82 

132.54 

183.20 

0.37 

1.33 

0.38 

0.65 

9.64 

7531. 

177. 75 

252.97 

134.47 

184.38 

0.40 

1.50 

0.33 

0.65 

9,64 

7475. 

179.02 

253.  12 

13  6.38 

185.  55 

0.  42 

1.67 

0. 33 

0.65 

9.64 

7418. 

180. 29 

253. 26 

1  38 . 29 

186.  72 

0.44 

1.83 

0.33 

0.  65 

9.64 

7362. 

181.54 

253. 41 

140.17 

187.88 

0.46 

2.00 

0.38 

0.  65 

9.64 

7305. 

182.  73 

253. 55 

142.05 

189.02 

0.43 

2.17 

0.  33 

0.65 

9.64 

7249  . 

184.  02 

253. 69 

143.90 

190. 16 

0.51 

2.33 

0.  33 

0.65 

9.64 

7193. 

185.24 

253.33 

145. 75 

191. 29 

0.53 

2.50 

0.33 

0.  65 

9.64 

7040. 

186.  45 

254. 21 

150. 67 

192.24 

0.  60  ■ 

2.67 

0.44 

0.  66 

6.  32 

3200. 

137.74 

251.15 

150. ei 

195.  60 

0 . 591 

2.83 

0.  50 

0.  67 

4.  03 

9645. 

169. 26 

246. 63 

146. 19 

200. 13 

0.53 

3.00 

0.  56 

0.  69 

2.51 

11241. 

191.03 

240.90 

143. 75 

205.80 

0.  S0 

3.17 

0.63 

0.  72 

1.52 

1  2373. 

193.07 

234. 08 

144.83 

212. 45 

0.  51 

3.33 

0.70 

0.  74 

0.9  1 

14368. 

195.  39 

226.93 

150. 73 

2 1 9 . 59 

0.  60  , 

3.50 

0.77 

0.78 

0.53 

15629. 

197.93 

220.47 

161.31 

226. 49 

0.77 

3.67 

0.85 

0.82 

0.31 

16501. 

200. 66 

215. 61 

1 74.88 

232. 49 

i  .ea  ■ 

3.83 

0.92 

0.86 

0.  13 

16992. 

203. 50 

212. 76 

189.14 

237.25 

1.53 

4.  00 

1.00 

0.90 

0.  10 

1 7201. 

206.40 

21 1 . 52 

201.93 

240.99 

2.  16 

4.17 

1.00 

0.90 

0.  10 

16922. 

209  .  £9 

213.  18 

203.94 

242. 77 

2.27 

4.67 

1 . 00 

0.90 

0.  10 

15647. 

217. 59 

220. 37 

212. 47 

246.21 

2.36 

CO  YOU  NEED  TO  STOPE  THE  OUTPUT  IN  A  DATA  FJLE?Y  OP  N 
srv 

EMTEP  A  5-LETTEP  FILE  NAME 
=  DLK  1 

YOMH  OUTPUT  VALVES  APE  STOPED  IN  FILE  DLK 1 

ENTEP  P.EPVN  OP  STOP 
=  STOP 
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REF:  ER  1110-1-10  -  ENGINEERING  AND  DESIGN  -  Engineering  and  Computer 
Program  Library  Standards  and  Documentation,  Appendix  C 

PART  III:  FILE  DOCUMENTATION 

1.  REVISION  LOG:  N/A 

2.  TITLE:  H5320  -  Lock  Filling  and  Emptying  -  Symmetrical  Systems 

3.  SOURCE  PROGRAM  LISTINGS:  See  pages  4l-52 

1*.  NUMERICAL  AND  LOGICAL  ANALYSIS:  A  finite  difference  scheme  is  used 
to  determine  the  flow  rate  through  the  culverts  and  water  surface  ele¬ 
vation  in  the  chamber  during  a  lock  operation. 

5 .  SUBROUTINES  NOT  DOCUMENTED  IN  ABSTRACT:  None 

6.  MISCELLANEOUS :  The  program  is  part  of  the  CORPS  computer  system. 
CORPS  is  an  acronym  standing  for  Conversationally  Oriented  Real-Time 
Program-Generating  System.  The  program  is  now  operational  on  the  WES 
G635,  Vicksburg,  MS;  HIS  66/80,  Macon,  GA,  and  Boeing  CDC,  Seattle,  WA. 
The  source  listing  on  page  4l  contains  the  first  line  run  command  and 
brief  for  H5320.  This  first  line  run  command  runs  the  binary  H5320B  of 
the  source  listing  on  pages  42-52  (Fortran  source  of  H5320)  and  attaches 


the  WESLIB  routines  HACCT  and  RERUN. 
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06/22/73 


0000U#RUN  *=J VESLIB/HS320B,  RCNOGO) 

10000  CHARACTER* 5  HFILE 
10010  COMMON  LOZ,LQX 
10020  HFI LEC SHH  5320 
10030  LQZ*l)LCX=l 
10040  PRINT  I® 

10050  10  FOPMATC /"H  5320  -  LOCK  FILLING  AND  EMPTYING- SYMMETPI  CAL  SYSTE 
10060  4MS"//> 

10070  15000  CONTINUE 
1003  0  CALL  HACCT< HFILE) 

10090  IFCLQZ.EQ.l)  CALL  H 53 20( $  1 6000) 

10100  GO  TO<20,  30),LOZ 
10110  30  CALL  DATAFC $16000) 

10120  CALL  MAIN 1< $16000) 

10130  PRINT,”  ”  1 

10140  20  LOZ=£ 

10150  CHARACTER  ZZZZZZ*2 

10160  16000  PRINT,  ”EN TE°  PEPUN  OP  STOP" 

10170  READ  16001,  ZZZZZZ 
10180  16001  F0FMATCA2) 

10190  I FCZZZZZZ . EO. 2HRE)  GO  TO  1500e 
10200  IFCZZZZZZ.E0.2HST)  GO  TO  20000 
10210  PRI NT, " EFFO R  ***  PETYPE" 

10220  GO  TO  16000 
10230  20000  STOP) END 

2000e  SUEP.OUTINE  H5320C*)  .  •*' 

20010  COMMON  LCZ  ,  LCX, XBBC  201),  TYPE,  R”N,XN,  Z  V,  ZL,  Z  P,XL  1  ,XL  2,  XL  3,  XL  4,  XL  5, 
200204KK 1,XK2,XK3,XK4,  DTS,  M,  I  VAL V, N,  BB 1 , AAB,  TV  1  ,  T 1  ,  TLAG,  T’.»2, AAC, SAG, I  PR  ' 
20030  4IN1,  IPP1N2,  I  DEL  1,1  DEL 2,  N  STAPT,  KK,  I VALl,  I  VAL  2,  LI  NEC  201  ) 

20040  DIMENSION  XTC  201),\KVC  201),XCC( 201), Q( 201 ), ppc £0 l ) , Z C  20 1 ) , Z UVC 201) 
200504, Z VC  201 >,Z DVC  201 ), v<  201 ), VHC  £01>,HLlC  201 >,HL2C  201) ,HLVC  201) ,HL3<  20 
200604  1) ,  HL4C 20 1 ),HLTC  £01),HI C  201 >,HC  201),XKC  201), Cl C  £01), VHVCC 201 ) 

20070  CHARACTER  RUN* 3 5,  AAA* 1 , AAB*8 , AAC* 1 0,  AAD*8, AAE* 1 0, CFILE+25 
20080  DATA  E, C, D, E, F, Cl C 1),XKV< 1 > /2. 2, 3. 2, . 65, . 3, . 9, 9 . 99, 1 0000.0/ 

20090  KK=0I  I  VAL  1  =  0)  I  VAL  2=0 

20100  10090  FOPMATC 3x, 3F3 . 3/3x, 5F3. 3/3X, AF8. 3/3X, F5. 1, 21 4, 5l 2) 

20110  10100  FORMATC 3x, 2F4, 0, A35) 

201£e  10110  FOPMATC 3x, 1 0F6. 3) 

20130*10120  FORMATC  "CON  DU  I T  GEOMETRY:  LENGTH  =  ”,F7.1,"  FT.  VI  DTH=  ", 

20  1  404  F  5.  1 ,  "  FT.",/,  13X,"HEIGHT=",F5.  1,"  FT.  N  =  ",l2, 

20  1  504"-CVLVEoT  OPEPATION)", //,  "CHAMBER  GEOMETRY:", 

201604-  LENGTH  =",F6.1,"  FT.  VI  DTH  =",F5.1,"  FT.  ", //,  "LO  SS  ", 
201704-COEFPIC1 ENT:  INTAKE", F7. 3, /, 9x, "UPSTREAM  COMDUI T", F7. 3, /, 

2018047k, "DOVNSTPEAM  CON DUI T", F7. 3, /, 1 7x, "MANI  FOL D”,  F7 . 3, /,  3x,  "TO TAL  . 
201904C  T/ALVE  OPEN)  ",  F7.  3,”  LOCK  COEF.  =”,  r  5.  3, //, 

20 2 00 4 "ELEVATIONS:  UPPER  POOL", F7. 2, /,  1 3x, "LOVEP  POOL", 

2021G4F7. 2,/,  1  lx, "CULVERT  POOF", F 7. 2) 

20220  10210  FORMATC //, "VALVE  OPENS  AT", F7. 2, "-MI N.  RATE  F0R",F7.2, 
202304"  MINUTES") 

20240  10230  FORMATC "VALVE  HOLDS  FOP", F6. 2, "-MIN.  AT", F6. 2, "FT. • ", /, 
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H5320S  CONT. 


20  2  50  4"  THEN  OPENS  AT’*,  F6.  2,  "-MI N .  PATE  TO  FULL  OPEN",/, 

202604 "TOTAL  VALVE  TIME  =", F6. 2, "-MI N.  C Nm", I  3, ">"//> 

20270  10260  FOPMAT</, "SPECIAL  VALVING,  TIME3",  F7.  2,  "-MIN.  "//) 

20280  10270  FORMATC  "EXTREME  ELEVATION  AT  TIME  «", 

202904F6.2,"  MI MUTES", /, "MAX  COMPUTED  OVERTRAVEL  *",FS.2,"  FT."> 

20300  10290  FORMATC//, 6x,  "VALVE  PARAMETERS", I 5x, "ELEVATIONS  < FT- DATUM) "/ 
2031 0S"TI ME  OPEN.  CONT.  LOSS  TOTAL  LOCK  US  VALVE  PIEZ.HD  DS", 

203204”  VALVE"/"  CM)  RATIO  COEF.  COEF.  INFLOW  CHAMBER  WELL" 

203304,"  VALVE  WELL”,/,  27x,"CFS",33x,"CAV. PAR.") 

20340  10330  FOPMATC F6.  2,  2F5.  2,  F9. 2, F7. 0, 4f8 . 2, FS. 2) 

2035e  10340  FORMATC  9x,  "RATE",  10X,  "VEL  HDSC  F)  ",  10X,  "DI  FFEPENTIAL  HDSCF)"/ 
203604, 9x, "PI  SE  VEL.  REF.  VENA  VALVE  TOTAL  INERTIA  OVERALL"/  ' 

203704, 8x, "C  F/M)  CF/S)  AREA  CONT. ",  1 lX, "LO SS"/) 

20388  10370  FOPMATC F6. 2,  3f7. 3, 3F8 . 3, F7. 3, f3 .3) 

20390  10380  FORMATC//, "LOCK  EMPTIES  TO  LOWER  POOL  IN  ",  F5.2, "  MINUTES") 
20400  10390  FORMATC//, "LOCK  FILLS  TO  UPPER  POOL  IN  ",F5.2,"  MINUTES") 
20410  10400  FORMATC "ACC. VEL.  ERROR  TO  I=", I  3, "  IS", F>0.5) 

20420  10410  FORMATC//) 

20430  10420  FORMATCAl, A8, Al> 

20440  10430  FORMATC  I  3,  3F6.  2,  F9. 2, F7. 0, 4F7. 2, 4F6. 2, Sf7. 2> 

20  450  N?TART=0 

20460  PRINT,  "DO  YOU  HAVE  A  PERMANENT  DATA  FILE?Y=YESi N=NO" 

20470  READ, AAA 

20430  IFCAAA.EO. lHY)GO  TO  10540 

20490  PRINT,  "WE  WILL  HELP  YOU  SET  UP  YOUR  DATA  FILE- 
20500  PRINT,  "YOU  MUST  SUPPLY  29  DATA  ITEMS  AND  A  FILE  NAME" 

20510  PRINT, "ENTER  YOUR  NEW  5-LETTER  FILE  NAME” 

20520  READ,  AAB 

20530  CALL  DATAFC $22275) 

20540  GO  TO  10610 

20550  10540  PPIMT, "ENTER  YOUR  5-LETTER  FILE  NAME" 

20560  N START3 1 

20570  10560  READ, AAE 

20575  IFCAAB.E0.5HST0P  )  RETURN  1 

20580  GO  TO  10610 

20590  10580  PRINT  10590, AAE 

20600  10590  FORMATC  "NO  FILE  ",  A0,".  ENTER  FILE  NAME  OP  STOP") 

20610  GO  TO  10560 

20620*10610  EMCO DEC AAC, 10420)  "/", AAB, "I ” 

2e63C  ENTFY  MAINlC*) 

20640  10630  CONTINUE 

20650  CALL  ATTACHC 1,  AAC,  3,0,1 STAT,) 

20660  CALL  TESTC I STAT, I  DEC) 

£0670  GO  TOC  10670,  10  530,941,941, 10530), I  DEC 

20680  10670  CONTI NUE»KK=1 

20690  PEA DC  1 ,  10102)  TYPE,XN, RUN 

20700  PEADC  1,  1  0090)  £*',  ZL,ZR,XLl,  XL  2,  XL  3,  XL  4,  XL  5, 

207  104XK  I,XK2,XK3,XK4,  DTS,M,N,  I  VALV,  I  PPI N  1 ,  I  PPI N2,  I  DEL  1 ,  I  DELS 
20720  I FC ABSC TYPE)  .LT. 0.99.OP. ADSCTYPE) • GT. 1 . 01 > GO  TO  294 
20730  IFCIVALV.EC.2)G0  TO  10750 
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H5320S  CONT. 


20740  READt l,  10110, END*10770) tXBBt I),I*1,N> 

20750  GO  TO  10770 

20760  10750  READt  1,  10110,  END- 10760)  SAG,  TVl,Tl,TLAG,TV2 
20770  10760  CALL  VAL2< SAG, TVl, T»,  TLAG, TV2, DTS,N,  EB»,XBB) 

20780  10770  IFtNSTART. EQ.0)GO  TO  10320 

20790  PRINT, "DO  YOU  WISH  TO  CHANGE  ANY  DATA  BEFORE  YOU  START?  YES  OP  NO’ 
203001" 

20810  READ, AAA'  I F< AAA, EQ« lHY)  KK» 1 
208  20  IFtAAA.EO. lHY)  CALL  RUN  1 C S 1 0630) 

208  30  103  20  CONTINUE 
20840  NS  TART5*  0 

20850  DTM=*DTS/60.  0J  t3=»(N-  1 )  *  DTMl  NC'XNi  AR=  CXN*XL2*XL3> ZtXL4*XLS) 

203  60  DO  109  30  1  =  2,  m+1 

20370  I  Ft  I .GT.N)XBBC I >=XBB< I- 1) 

20330  X  Tf  I  )  “X  Tt  I  -  1 )  +  DTM 

20390  I FCXBBt I ) .LE. 0. 2)  XKVt I > =t 1 0. 0** t B- . 2*C> ) * t 0. 04/XBBC I ) **2) 

20900  I FCXBBt I ) .GE.0. 2)  XKV< I >“10.0**t B-C*XBBt I > ) 

20910  IF<XBECI).LE.0.3>  XCCt I ) - 1>  t E- D) *  CO St  3.  1 4 1 59 3* (XBB(  I ) /0. 60) ) 

20920  IFCXEE(I) .GE.0.3)  XCCt I )  =  F- 1 F-D) *CO St t  tXBBt  I  > - .  3)  * 3.  1 4 1 59 3) / 1 . 4> 
20930  XKt I ) =XK 1+XK2+XK3+XK4>XKV< I ) 

209  40  109  30  CONTINUE  !  ‘ 

20950  XKT=XK1+XK2+XK3+XK4*10.0**<B-C) 

20960  XKC=SOPT( 1/XXT)J PRINT  10950, RUN! 10950  FOPMATC /a35//) 

20970  IFClPRINl.EQ.  DPRINT  1  0  1  20, XL  1,XL2,xl3,NC,XL4,xlS, 

20980«XKl,XK2,XK3,XK4,XKT,XKC,ZU,ZL,ZR 

2099 C  BB2=XL3*EBl 

21000  I  Ft  I VAL V. EG.  2) RRI NT  10210,  TVl,T> 

21010  IFtIVALV.E0.2.AND.Tl.LT.TVl)PRINT  1  0230,  TLAG,  EB2,  TV2,  T3,N 

21020  I  Ft  I VALV. EG. 1 >PPINT  10260, T3 

2103ZC  SET  INITIAL  CONDITIONS 

21040  J Ft  TYPE) 1  1040,  1 1060,  1 1 060 

21050  1  1040  Et  1  )=£!»;  ZPEF*ZLiZDWt  1)*ZU 

2ie60  GO  TO  11070 

21070  1  1060  Zt  l>=ZL«'ZPEF*‘ZU'ZDVt  l)=ZL 
21030  11070  CONTINUE 

21090  XCCf  l)“ElMl=MlM2*MlZUWt  D«ZUiZOT*0.0IZVt  l)mZL 
21 100  Mo=M 

21110  Ht  1 )  =ZU-ZL»  HL  Vt  1  >**Ht  1  )  IHLTt  I  >  *Ht  I  ) 

21120(?  STAPT  NUMERICAL  INTEGRATION 
21130  DO  11270  I»1,M-1 

21140  Vt  1  + D=r*’EDf  Z  f  I  ),  Vt  I  >,XKf  I  ♦  1  ),  TYPE,  ZREF,  AP,  DTS,XLl) 

21150  Zt I ♦ 1 ) “ZNEVt  Z  t I ), Vt I ) , Vt !♦! ), AP, DTS, TYPE) 

21  160  I  Ft ABStXL 1* t  Vtl  ♦ 1 ) - Vt I ) ) > .  GE. . 5*XKt  I  ♦ I ) * DTS* Vt  I  ♦ 1 > ** 2>m3  =  I -  1 
21170  IFf I.GT.M3.0P.I.GT.Ml)G0  TO  11250 
21130  CO  11230  J«l,  10 

21  190  Vt I ♦ 2) =PPEDt  Z<I  +  I),V(I*1> ,XKt 1*2), TYPE,  Z  PEF, AR, DTS, XL  1) 

21  200  v< I ♦ 1 ) -COPRf  Z 1 1 ) , Vt I ) ,XXt I ♦ 1 > , Vt I ♦ 2)  ,  TYPE,  Z PEF,  AP,  DTS, XL  1 ) 

21210  ZCHEC*ZNEVf  Z t I), VtI),Vt  !•*■!),  AP, DTS,  TYPE) 

21220  IFtABStZCHEC-Zt !♦!)). LE.0.0001)GO  TO  11240 
21230  Z 1 1 ♦ 1 ) =Z  CHEC 
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21240  11230  CONTINUE 
21250  11240  ZC I ♦ 1 )«ZCHEC 

21260  11250  IFCM2.E<5.M.AND.V<I+1>.LT.0.0>M2«I 

21270  lFCMl.EQ.M.ANC.TYPE*<ZREF-ZCI«-l>>.LE.0.6>  Ml*I 

21280  11270  CONTINUE 

21290  EPPOR=0.0 

21300  PRINT  10410 

21310  DO  11350  I  =  2, m2 

21320  YY=COPPCZCI-l>,  VC  I  -  »  >  ,  XKC  I  >  ,  VC  I  +  1  >  ,  TYPE,Z  REF,  AR,  DTS..XLI  > 

21330  EPROF=CYY-VC I > >**2+ERP0P 
21340  I FC I . EQ.M 1 >  PRINT  1 0 40 0, I , ERFO R 
21350  I FC I . EQ.M2) ppiNT  10400, I , ERPOR 
21360  11350  CONTINUE 
21  370C  END  OF  DI  FF •  EON  S* 

21330  DO  11690  l  =  2,M  1 

21390  VHC I ) =C VC  I ) ** 2} / 64. 4 

21400  OC I >  =  VC I )*XL2*XL3*XN 

21410  FRCI)=CeCI)*60.0)/<XL4*XL5> 

21420  HI C I)=XLl*C VC  I* l >-VC I  -  1) )/C  64. 4*DTS> 

21430  HLlCl)=XKUVHCI> 

21440  HL2(I)=XK2*VHCI>  : 

21450  HLVCI)=XKVCI)*VHCI> 

21460  HL3<  I  )=XK3*»7HC  I) 

21470  HL4c  I  )=XK4*tjHC  I  ) 

21430  HLTC  I  )  =XKC  I  )  *’fll  Cl) 

21490  HC I ) =TYPE*CZFEF-ZC I ) ) 

21500  IFCVCI)  .  6E.  0.0)VHVCCI>=VHCI>/CXCCCI>*XBBCI>>**2 
21510  IFCVCI) .LT- 0- 0) VHVCCI)=VHC I) 

21520  I  FC  TYPE)  1  1590,  1  1590,  1  1520 

21530  1  1520  IFCVCI)  .GT.0.0)ZWC1)=ZU-C  1.0+XKl+XK2)*VHCI) 

21540  IFC  VC  I  ) -LE.  0- 0)ZUVC  I  )=ZU-C  1  •  0-XK»-XK2)*  WCI  ) 

21550  IFC  VC  I  ) • GE. 0.0)ZDVC I ) =ZC I >+XK4* VHC I ) 

21560  IFCVCI) .LT.0.®)ZDVCI)=ZCI)-XK4*VHCI) 

21570  I FC VC  I > .GT. 0. 0>Z  VCI >=ZUWC I )+VHC I > -VH VCC I ) 

21  580  IFCVCI)  .LE.  0.  ®>ZVC  I  )=ZUVC  I  >  -  VHVCC  I  ) 

21590  GO  TO  11640 
21600  11590  ZUVCI)=ZU 

21610  I  FC  VC  I  )  .GE.  0.  01ZDVC  I  >=ZC  I  )  -C  1 . 0+XK  1+XK2)*  CMC  I) 

21620  «:  FC  VC  I  )  .LT.  £>•  e>ZDUC  I  >=ZC I > -C  1.0-XKl  -XK2)*VHC  I) 

21630  IFC  VC  I  )  .GE.  0.  0)ZVCI  >=ZDVC  I  >  +  VH  C  I  >  -  VHVCC  I  > 

21640  IFCVCI). LT.0.6)ZVCI)=ZI>UCI)-VHVCCI) 

21650  11640  CONTINUE 

21660  IFC VC  I > .LE. 0. 0)CIC I >=9.99 

2167  0  I FC VC  I ) . GT.  0. 0) Cl C I  )  =  C  Z VC  I  > -ZP+XL3*C 1 . -XCCC I )  +  XBBC I ) >> /VHVCC I ) 
21  680i-f33. 0/VHyCC  I  ) 

2169  0  I FC  Cl C 1 ) . GE. 1 0. 0) Cl C I > =9 . 99 
21700  11690  CONTINUE 

21710  XTX=XTCMl)+DTM*CZPEF-ZCMl>)/CZCMl+l>-ZCMl>) 

21720  XTM=XTCM2) 

21730  Z0T=ZCM2) -ZREF 
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21740  I  FCM  1  •  EQ.M) PRIN  T*  "LOCK  OPERATION  IS  NOT  COMPLETE” 

21750  IFCM2.E0.M)  PRINTS  "LOCK  DOES  NOT  COMPLETE  OVERTPAVEL" 

21760  l  FCMl  .LT.M.  AND. TVPE.GT. 0. 0) PRINT  10390.XTX 
21770  IFCMl. LT.M. AND.  TYPE. LT.0.0>PRINT  10380.XTX 
21780  I FC M2. LT.M) PRINT  1 0270.XTM* ZOT 
21790  PRINT  10290 

21800  PRINT  10330, CXTCI)*XEBCI),XCCCI  >  *XKVC  I)*Q(I)*ZCI)*Z  UVC  I  >* 
218I04tV<I)jZDW<I)^CI<I)/I*l>N,IDELl) 

218  20  PRINT  10330,CXTCI>*XBBC  I ) * XCCC I > *XK VC  I > ,  QC  I  >  * Z C  I  >  , ZUVC  I  >  * 

218  304ZVCn.ZDVCI)*CICI)*I=N+lDEL2,M,IDEL2> 

21840  IF  CIRRIN2)  1  1900.  11900.  1  1340 

21350  1  1840  CONTINUE 
21863  PRINT  10340 

213  70  PPINT  10370, CXT<I>jRP(  I  >  *  VC  I  >  *  VHC  I  )  ,  VH  VCC  I  )  ,HL  VC  I  >  *HLTC  I  >  , 
218804HICI)*HCI)*I=l*N*IDELl> 

21890  PRINT  10370,  CXTC  I  >*RRC  I  )*  VC  I  >*  VHC  I  >*  VHVCC  I  >*HLVC  I  )*HLTC  I  >* 
21900 4HICI).H(I),I=N+I DEL2.M, I  DEL 2) 

21910  11900  CONTINUE 

21920  PRINT,  "DO  YOU  NEED  TO  STORE  THE  OUTPUT  IN  A  DATA  FILE?Y  OP  N" 

219  30  PEAD,  AAA 

21940  IFCAAA.EQ. lHN)GO  TO  320 

21950  PRINT*  "ENTEP.  A  S-LETTER  FILE  NAME” 

21960  PEAD*  AAD 

21970  LI  NEC  1 )  =  1 00 

21930  DO  11980  I  =  2,  M 

21990  11980  LINECI)=LINECI-1)>2 

22000  ENCODEC  AAE*  1 0420)  "/•’*  AAD*  "S  " 

22010  ENCODEC  CFILE*  59)  "CF* /",  AAD,  ",  B/ l*  35/,  R,  W#" 

22020  CALL  ATTACH  C  2, AAE, 3, 0, I STAT* ) 

22030  CALL  TESTC  ISTAT*  I  DEO 

22040  CO  TOC  1 72,  1 73,  1 75,  1 75,  1 75) , J DEC 

22050  173  CALL  ACCESSC CFILE* $1 75) 

22060  172  CONTINUE 

22070  CALL  ATTACHC 2, AAE* 3, , ) 

22080  CALL  HEADERCM) 

22090  DO  12075  1  =  1, Mi  12075  XKVC I ) *ALOGl 0CXKVC I > ) 

22100  VPI  TEC  2,  1  0430)  C  LI  NEC  I  )  *XTC  I  >  ,XEBC  I  >  *XCCC  I  )  *XKVC  I  )  ,  CC  I  )  *  ZC  I  >  * 
22110*7 UUC  I  >  *  Z  VC  I  )  * Z  DWC  I  )  ,  Cl  C  I  )  *  RPC  I  )  ,  VC  I  >  *  'XU  I  >  *  VH  VCC  I  >  *HL  VC  I  >  , 
2212C««:lTCI)*HICI),HCI)*I',1*M> 

22130  61  CALL  DETACH  C  2* , ) 

22140  PRINT  160, AAD 

22150  160  FOPMATC "YOUR  OUTPUT  VALUES  ARE  STORED  IN  FILE  "*A8> 

22160  59  F0RMATCA4, AS, Al3) 

22170  GO  TO  320 

22180  175  PPINT*  "FILE  PROBLEM;  ILLEGAL  CHARlBUSYiOP  SYSTEM  PROBLEM" 
22190  PPINT* "OUTPUT  VALUES  NOT  SAVED  AS  REQUESTED" 

22200  320  CONTINUE 
22210  GO  TO  293 

22220  294  PRINT*  "PROGRAMMED  STOP*’ TYPE  NOT  EQUAL  1  OR  -1"*’PETURN  1 
22230  941  PRINT  946, ISTAT 
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22240  946  FORMATC "FILE  PROBLEM.  CALL  FTS  542-2403  AND  GIVE  THIS  NIMB 
222504EP  ",0l2> 

22260  29  3  CONTINUE 
22270  CALL  DETACH < 1 ,  ,  )  '  RETURN 
22275  22275  RETURN  1 
22230  END 

22290  FUNCTION  PPED  C Z  ST,  VST.XKF,  TYPE,.  Z REF,  AR,  DTS.XL I ) 

22300  XC  =  <  64. 4/XKF) *C  TYPE*  CZREF-Z  ST> - VST*C . 5*AP*  DTS-XL l/C 32. 2*DTS) ) ) 
22310  XB*C  64. 4/xKF) *< . 5*AR* DTS+XL 1/C 32. 2*DTS) ) 

22320  IFCXC-GE.  0. 0) PRED= 5*XB+ . 5* SQRTCXB** 2+ 4.*XC) 

22330  IFCXC.LT. 0. 0) PRED* . 5*XB- . 5* SQRTCXB** 2-4. *XC) 

22340  RETTJPM 1  END 

22350  FUNCTION  CORR  C Z ST, VST, XKF, VAH, TYPE.Z REF, AP, DTS.XL 1 ) 

22360  XC=C 64. 4/XKF) *C TYPE* C Z REF-Z ST) - . S* VST*AR*DTS-XL 1* C VAH- VST)/ 
223704C 64.4* DTS>> 

2233  0  XB*  32. 2*AR* DTS/XKF 

22390  IFCXC.GE. 0. 0) COPP=- . 5*XB+ . 5* SQRTCXB** 2* 4. *XC> 

22400  I  FCXC.LT.0.0) COPR=. S*XB-. 5*S0RTCXB**2-4.*XC» 

22410  RETURN) END 

22420  FUNCTION  ZNEV  < Z ST. VST. VF, AR, DTS. TYPE) 

22430  ZNEV=ZST+TYPE*.  5*AR*DTS*C  VST+VF) 

22440  return; END 

22450  SUBROUTINE  TESTC I STAT, I  DEC) 

22460  DIMENSION  I  TESTC  5) 

22470  DATA  I  TEST/000. 005.037, 004,034/ 

22430  I  DEC* 3 

22490  I  STAT=FLDC  6,  6,  I  STAT) 

22500  I FC  I  STAT. EG. I TESTC 3). OR. I STAT. EQ. I TESTC l>>  I  DEC* 1 

22510  I FC I STAT. EG. I TESTC 2) )  IDEC=2 

22520  I FC I STAT. EQ. I TESTC 4))  IDEC*4 

22530  IFCISTAT. E0.ITESTC5))  IDEC=5 

22540  RETURN  .’END 

22550  SUBROUTINE  VAL 2C A, TVl , Tl . TLAG. TV2, DTS.N,  BBl.XBB) 

22560  DIMENSION  XTC  20  1 )  ,XEBC  20  1  > 

22570  t2s=0.0;bb1*0.  0.‘  DIM*  DTS/6B.  0)XTC  l  >  =  0.  0)  PI*  3.  I  41 59  3 
22580  I  FC  TV  I  -  DTM)  132,  136,  186 

22590  13  2  PRINT, "YOUR  1  ST  VALVE  RATE  IS  LESS  THAN  ONE  TIME  STEP" 

22600  PRINT, "ENTEP  VALVE  OPENING  PATIO  AT  END  OF  I  ST  TIME  STEP" 

22610  Vf.AD,  BB 1 

22620  Tl*0.0 

22630  IS 6  CONTINUE 

22640  lFCTVl.GE.DTM)EBlaTl/TVl-A*SINCPI*Tl/TVl> 

22650  T2s=BBl*TV2 
22660  DO  214  J- l, 6 

22670  t2n*T2S*C  TV2*EE * -T2$+ A* TV2* S1NC FI *T2S/TV2>  > / 

226304C 1.0-A*PI*COSCFI*T2S/TV2> ) 

22690  I FC AESC T2N/T2S- 1 • 0) .LE. 0.0001)  GO  TO  2*5 
22780  214  T2S=T2N 
22710  215  CONTINUE 
22720  T2*Tl+TLAG 
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22730  t3*t2-t2c  +  t,’2 
22740  N  =  t3/ CTM+  £ 

22750  DO  52  1  =  2, N 
22760  XTC  I)=XTC  I  -  D+CTM 

22770  IFCXTC  I  )  .LE.Tl  )XEEC  I  )=XT<  I  ) /T”1-A*FINC  FI*XTC  I  )/TVl  > 

2273  0  I  FtvT( I ) . GT. T 1 . AND.  XTC I  )  .LE. T£)XEEC I ) =EB 1 

22790  I  FC  XTC I ) . GT. T2>XEE< I  )  •  C XTC I ) - t2+ t£? ) / TV2- A*  T i NC El  * CXTC I > -t2+T2F> 
22300*/T,,2) 

223X0  52  CON"rIN"E 

223  20  I F CXEEC  N) . GT.  1 . 0)XEEC  N  )  =  1 . 0 

223  30  rET,,r,N'END 

22340  e»*ET>0,,‘’INE  DATAF C  *  ) 

223  50  COMMON  LOT,LOX»XE£C  201  )  ,  T'PE,  E»'N,  XN,  T  »r,  Z  L,  Z  P,XL  1 ,  ML 2, XL 3,  vl4,  XLS, 
223  604XK 1 ,  VK2,XK3,  XK4,  CTe,  M,  I  VALV,N,  EE  1, AAE,  TVl,  Tl,  TLAG,  TV2, AAC,  FAG,  I  PP 
£23  704INl,Ic'cIN2,I  DEL  1 ,  I  EEL2 ,  N  CTAPT»  XK,  I  VAL  1 ,  I  ”AL2,L1NEC  2&1>  ' 

2238  C  CHApAC'rE°  diN*  35,  AAE«3  ,  AAC*  1  0,  CFI  LE*  2  A 
22390  I FCLOX. EO. 3)  GO  'TO  X 1 520 
22900  £>r™  «>"Nl(*) 

22910  I FCLCT . EO. 2)  GO  TO  11520 
22920  CALL  PPINTT 

£2930  I  F C *J CTAT?T«  EO .  0 )  GO  TO  1  1520  ,  •’ 

£2940  LOZ  -2*  PEVIND  1 

22950  I  1520  KKK  =  29 

22960  GO  TOC  1  I  1,  *  12), LOT 

£2970  1  12  CALL  PE^’NC  KXK,  LCX,  JKL  ) 

2293  0  GO  TOC  1,2, 3, 4, 5, 6, 7, 3, 9,  10,  11,  12,  13,  14,  15,  16,  17,  18,  19,20,21,22/23, 
22990*24,  25,  26,  27, 23, 29, 30) , JKL 
23000  111  CONTINUE 
23010  GO  TOC  1  ,  1  1  2)  ,  LOT 

23020  1  '■'’I'JT,  "AA-EN~Eo  1=FILL^0P  -  1  =  EMP'T''  OP  0  =  TECMI  WflTE  NOV” 

23030  PEAD,  •'~VT'E 

23040  I FC APFC TY”E) .LT. 0.99) GO  TO  174 
23050  GO  TO C 2,  112), LOT 

23060  2  °PI!1"‘,  "AE-EN7E0  Y0”0  PPODLEM  NAME” 

23070  "EAD  33,  P”N 
23030  33  FOPMATC  A35) 

2309  0  GO  'rOC  3,  11  2)  ,  LOT 

23100  3  ppIN~,,,AC-EN-rEP  1  =  1  -CPLV.  'OP  £=2-C'’LV.  OPEPATION" 

23110  EAD,  XN 

23120  GO  tq< 4, 1 1 2) , LOT 

23130  4  POINT/  "AD-EM-t-e^  ,,t’P EP  FOOL  ELEVATION” 

23140  PEAC,’*’ 

£3150  GO  TOC  5,  1  12), LOT 

£3160  5  PPI NT, "AE- EN TED  LOVEP  POOL  ELEVATION” 

23170  "EAC,TL 

£3I3e  GO  T0C6, 1 12), LOT 

23190  6  nni«r, "AF-EMTEF  CVL".  POOF  ELEVATION" 

£3200  PEAD, TP 

23210  GO  TO C7,  112), LOT 

2322G  7  P"lMf,  "AG-ENTEP  C"L".  LENGTH” 
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23230 

23240 

232S0 

23260 

23270 

2328  0 

2329  0 
23300 
23310 
23320 
23330 
23340 
23350 
23360 
23370 

2338  0 

2339  0 
23400 
23410 
23420 
23430 
23440 
23450 
23460 
23470 
2343  0 
2349  0 
23500 
23510 
23520 
23530 
23540 
23550 
23560 
23570 
23580 
2359  e 
23600 
23610 
23620 
23630 
23640 
23650 
23660 
23670 
23630 
23690 
23700 
23710 
23720 


PEAD,  XL  1 

GO  "Oca,  1  12), L gz 
3  PPINT,  "AH-ENTEP  C*’LV.  %’IDTH” 

READ/XL2 
GO  ‘r0<9/  1 1  2), LOT 

9  PPINT,  ”AI  -ENTEP  C’’LV.  HE!  GTH” 

PEAP,VL3 

GO  TOC  10/ 1 12), LOT 

10  PPI  NT/  ”AJ -  EN  TER 
°EAD,  XL4 

GO  to<  1  I,  1 12), LOZ 

11  PRINT,  "AK-ENTER 
PEAD,XL5 

GO  TOC  12/ 112) /LOT 

12  PPINT,  ••AL-ENtE'° 

PEAD/VKl 

GO  TOC  13, 1 12), LOT 

13  PPINT, "AM-ENTER 

READ/XK2 

GO  TOC  14,  11 2), LOT 

14  PPINT, "AN- ENTEo 
PEAD/VR3 

GO  TOC  15,  1  12), LOT 

15  PRINT, -AO-ENTEP 
PEAD,  XX  4 

GO  toc 16, 112), LOZ 

16  PPINT,  "Ar- ENTEP 
PEAD,  DTP 

GO  TOC  17, 1 12), LOT 

17  POINT, "A0-ENTE0 
PEAC/M 

GO  toc 18,  1 1 2)  , LOZ 
13  PPINT,  •'A0- ENTEP 
PEAD,  l"AL" 

IFCIVAL”. E0.2)  GO  TO  50 
GO  toc  19,  1  1556), LQZ 
11556  LOZ=l 

19  PRINT, "AS-ENTEP  NO.  OF  PATIO?  TO  EE  ENTEPED" 
•PEAD/N 

GO  TOC  20, 1 1 2) ,L0Z 

20  ppint, "AT- ENTEo  paLVE  PATIOSCSEP.  EY  COMMAS)” 
PEAD,  CXBBC  I  >,  I3  l,N> 

I FC I VAL 1 • EC • 0. AND. 1 VAL2. EO • 1)  LOZ«2 
I  *'AL  1 3  1 

IFCKK.E0.1)  LOZ 3  2 
GO  TO  60 

50  GO  T0C21, 11591), LOZ 
11591  LOZ3 1 

21  ’’’’INT,  "A”- ENTEP  ”ALVE  SAG" 

PEAD,  SAG 


CHAMBER  LENGTH" 

CHAMBER  WIDTH” 

CI’LV.  INTAKE  LOSS  COEF.” 

v-s  cav.  LOSS  COEF.” 

D-S  CI’LV.  LOSS  COEF." 

C”LV.  EXIT  LOSS  COEF.” 

TIME  STEP  IN  SECONDS" 

TOTAL  NO.  OF  CALC.  STE”?  DESIFED” 
TYPE  OF  VALVINGC  1  OP  2)  FEE  TABLE) 
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23730  GO  TO<  22,  1  1  2)  , LOZ 

237*0  22  PPINT,  "AV- ENTE®  FI  ®ST  VALVE  PATE" 

23750  PEAD, TV I 

23760  GO  TOC23, 1 12),L0Z 

23770  23  PPINT, "AV-ENTE®  TIME  VALVE  STOPS  MOTION!  MINVTES) " 

2378  0  PEAD,  Tl 

23790  GO  TOC  24, 1 12) ,LQZ 

23800  24  POINT,  "A'C-ENTER  DVPAT1 ON  OF  STO  P<  MI  NVTES)  " 

23810  PEAE,  TLAG 

233  20  GO  TOC  25,  112), LOT 

23330  25  POINT,  "AY- ENTER  FINAL  VALVE  PATE" 

233  40  PEAC,  TV2 

23350  IFCIVALl.EO.  1.AND.IVAL2.EO.0)  L0Za2 

233  60  I'»al2*1 

23370  IFCKK.E0.1)  LCZ®2 

23330  60  CONTINUE 

23390  GO  TOC  26,  112)  ,L<?* 

23900  26  poiNT,  "AZ-ENTEP  1ST  Ot’TPVT  CONTPOL*’  l* DO,  0*  CO  NOT  PPINT" 
23910  ®EAC, ioojn  1 
23920  GO  TOC  27, 1 12) ,LQZ 

23930  27  PPINT, -BA-ENTEo  2ND  OVTPVT  CONTPOL' 1- CO, 0“ DO  NOT  PPINT" 
23940  PEAE,  Iooin2 
23950  GO  TOC  23,  1 12),LQZ 

23960  23  PPINT, -be- ENTEP  El  pop  Ot»TP"T  PATE'*  1®E’'EPY  VALVE" 

23970  PEAE, I  CEL  1 
23930  GO  T0(29, 1 12),L0Z 

2399  0  29  PPINT,  "EC-ENTEP  FINAL  OVTP’.’T  PATE>2*EV.  2ND  VALVE" 

24000  ocad,  I  DEL2 

24010  GO  TOC  30,  1 1 2) ,LCZ 

24020  30  CONTINUE 

24030  LI  NEC  1 ) = I  00 

24040  CO  54  J=2,?J/10*6 

24050  54  LINECJ)=LINECJ-l)+2 

24060  ENCOCECAAC,  55)"/", AAE, "7 " 

24070  CO  62  J=M+1,N+10 
24030  62  XEEC J) "XEBC J- 1 > 

2409  0  ENCODE!  CFILE, 59>  "CF, /", AAB, ”, B/l,  3/,P,V#" 

24100  CALL  AttacNC 1 , AAC, 3, 0, I STAT, > 

241 ie  CALL  tEST( I STAT*  1  DEC) 

24120  GO  TOC  1  72,  1  73,  1  75,  1  75,  l  75),  I  DEC 
24130  173  CALL  ACCESS! CFI LE, I 1 75) 

24140  172  CONTIN"E 

24150  CALL  ATTACH < 1 , AAC, 3, , ) 

24160  t.rpi  TEC  1 , 5  6)  LI  NEC  1  )  ,  TYPE,XN,  PVN,LINEC  2>,^f,  ZL,ZP,LINEC3>, 
241 7eiXLl,XL2,XL3,XL4,XL5,LINEC  4) ,XK 1 ,XK2, XK3,XK4, LI  NEC  5) , DTS, 
2413  0<M,  N,  It'AL”,  IPPIN  I,  IP®IN2,  l  DEL  1 ,  l  DEL  2 
24190  DO  53  J=1,N/10+1 

24200  I  FC  I’fALV.  EQ  ♦  2 )  v®l  TEC  1 , 5  7  )  L I  NEC  6)  ,  SAG,  TV  l ,  T 1 ,  TLAG,  TV2 
24210  iFCI’fALV.  E0.2)C0  TO  61 

24220  58  opItec 1 , 57 )L IN  EC J* 5) ,CXBEC I >, I* 10*3-9,  10*J> 
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24230  61  CALL  DETACH ( Ijj ) 

24240  55  F0RMAT!Al.A8,Al) 

24250  59  FOPMAT! A4, A3,Al2) 

24260  56  FORMAT!  t  3,  2f4.  0,a35,/,I  3,  3f8. 3,/,I  3,  5F8. 3,/,  I  3,  4f3. 3, /, 
24270*1 3, F5. 1, 2l 4, 51 2) 

24230  57  F0PMAT!l3, 10f6.3> 

24290  GO  TO  200 

24300  174  PPINT,  "PPOGPAMMED  STOP' tyPEk0»  PEVI EV  VO”?  DATA"!  PETUPN  » 
24310  175  oointfILE  PPOELEMS'"  ILLEGAL  CH APACTEF'  BUSY'  0 F  SYSTEM  PPOB." 
24320  200  IFCNSTAPT.EO. 1)  GO  TO  l 192 
24330  PET?'PN 

24340  1192  N  START* 0' L02  *  1 ' PETUPN  1* END 
24350  S»'BP01'TINE  »PINTT 
24360  CHAPACTEP  AAA* I 

24370  PPINT,  "SHOfT-D  I  PPINT  A  TABLE  OF  INPUT  CO DES?Y=YES'  NaNO” 

24330  PEA D, AAA 

24390  IF!AAA. EO. lHN)GO  TO  4 

24400  PPINT  2 

24410  2  FORMAT! "CO  DE  VAPI  ABLE", 1 7x, "DESCPI PTI ON"/ 

244204"AA  TYPE’S  3x,’’lsFILL'-  1“ EMPTY'  0“  TEPMI MATE  P*TJ”/ 

244304”AE  PUN", 4x, "PROBLEM  NAME!  35  epACES  AVAI LAELE) ”/ 

244404’’AC  XN’S  5x.  *’l=ONE  C'*LV.J2«TV0  C'tVEPT  OPEPATION"/ 

244504-AD  Z*”S  5X, "UPPEP  POOL  ELEV. ! FT- DAT”M) ”/ 

244604-AE  ’L", Sx , "LOVER  POOL  ELEV. ! FT-DATVM) ”/ 

244704-AF  ZF",  5x,"C,tLV.  POOF  ELEV.  <  FT-  DAT* W>  AT  VALVE*’/ 

244304-AG  XLl’S  4x, "CVLV.  LENGTH<  FT) ”/ 

244904-AH  XL2'S  4y,  "CULV.  VI DTH< FT) "/ 

245004-AI  XL3’S4X,”CULV.  HEI  GWC  FT)  "/ 

245 1  04’’AJ  XL4",4x,"CHAMEEF  LENGTH!  FT)  ”/ 

245204) 

24530  pplNT  3 
24540  3  FOPMAT! 

245504-AH  YL5'S  4x,  "CHAMEEP  VIDnUFT)”/ 

245604-AL  XK 1 ”, 4X, ”1 NTAK E  LOSS  COEF.’V 
2457 04 ’’AM  XX2’S  4x,”»’-S  COND'UT  LOSS  COEF.'V 
2453 04 ’’AN  XK3’S  4X.  "C-S  CONDUIT  LOSS  COEF.'V 
245904-AO  XK4", 4x, ”EXI ?  LOSS  COEF.'V 

246004/ ’'NOTE! LOCK  COEF.=SO”APE  POOT  OF  l/CXK  1+XX2+XK3+XK4+0.  1 )  "// 
246104-AP  DTS”, 4X, "TIME  STEP! SEC.)  FOP  CALC”LATI  ONS’V 
246204-AO  M",6x,”N0.  OF  STEPPED  VALUES!  I NCL .  TIME=*0.0)’V 

246304-AP  IVALV’S2x,”VALVE  OPTI  ON  S:  1  =  PO I  NT- BY-POINT”/ 

246404.  25x,  ”2«U?E  PPOGPAMMED  FUNCTIONS’’/ 

246504/-IF  IVALV=>”// 

246604-AS  N’S  6x. -NO.OF  VALUES  OF  VALVE  OPENING  PATIOS  TO  EE  PEAD’V 
246704,  l  lx,’’!  INCL.  0.0  AT  TIME=0.0  AND  1.0  AT  VALVE  FULL  OPEN)’’/ 

2463  04,  1  lx, ’’THE ce  VALUES  APE  DTS  SEC.  APART!  SEE  AP,AEOVE)’V 
246904-ax  XBE! I ) ",  IX, ”N  VALUES! SEPAPATE  VI TH  COMMAS)  OF” 

247004/,  1  lX,  "VALVE  OPENING  PATIOS’’/ 

247 1 04/-I F  I  VALV=2”// 

247204-AV  SAG’S  4X, "NO PMAL  SAG  AT  1/2  VALVE  TIME!  APPROX. “0. 1 > ”/ 
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247304) 

24740  PPINT  6 
24750  6  FOPMATC 

247604-A"  TVl",4x, "INITIAL  OPENING  PATECMINS.  TO  FULLY  OPEN)"/ 

247704"A»*  Tl”,  5x, "VALVE  MOTION  STOP?  AT  Ti  MINE."/ 

247804"AX  TLAG",  3X,  "VALVE  HOLD?  P0  5ITI0N  FOP  TLAG  MIN?.”/ 

2479 04"AY  TV2", 4X, "FI NAL  OPENING  PATECMINS.  TO  FULLY  OPEN)"/ 

243 004 /"NOTES  NON- STEPPED  OPENING'  Tl* TVt  ,  TLAG* ®.  8,  TV2»TVl"/ 

243  1  04/"NOTEs  THE  NEXT  ITEMS  CONCEPN  THE  PPINTED  O^TPVT"/ 

248 204" At  IPPINl”, lX,"l*PPINTS  OUT  ITEMS  AA  THPV  AO, ABOVE”/ 

243304,  1 lx,  "0=DO  NOT  PPINT  ITEMS  AA  THPV  AO"/ 

243  404”BA  IPPIN2", lx, "LOSS  Cl  STPI  BUTI ON?i  1»PPI NT,  0*DO  NOT  PPINT"/ 
243504"EB  I  DELI  PPINT  1ST  AND  EV.  I DELl*DT?  EEC  TO  VALVE  F"LL  OPEN"/ 
243604"EC  ICEL2", 2X,  "PPINT  EV.  I  DEL £* DTE  SEC.  TO  <M-1)*DTS  EEC.”//) 

243  7  0  4  CONTINME 
2438  0  PET»'°N»END 
2489  e  Sf’EPO  VTI NE  HEADEPCM) 

24900  14840  FOPMATC ”08  H5320  18  15  EDGE"/"09  C 1 5< /) , C 3x, 3f6. 2, F9 . 2, F7. e, 
249  104 "77. 2,  4Fb.  2,  5f7.  2)  )  "/"10",  18  <  tX,  I  3) /”  1  l ",  18  <  ”  1”)  ) 

24920  14860  F0PMATC"12  C"°VE  DESIGNATIONS  FOP  H5320  APES"/"13  1*TIME”,2 

249 3044X," 10= CAVITATION  PAPAMETEP”/ " » 4  2*VALVE  OPENING  PATIO",  9X, ” 1 1 *PA 
249404TE-OF-PI  ?E"/"15  3=CONTPACT10N  COEF. ",  1 IX, ” 1 2« VELOCI TY  AT  PEF  APEA 
249  504"/" I  6  4=L0G<VAL"E  LOSS  COEF- >", 7X, " 1 3- VELOCI TY  HD  AT  PEF  APEA"/"1 
2496047  5*  TO  T  AL  1NFLOV",  1  t>X,  "1  4*VELO  CITY  HD  AT  VENA  CONT"/"13  6-LOCK  V 

249704ATEP  SVRFAC  EC  Y- S)  ”,  5x,  ”  1 5*H  EACLO  SS  THPV  VALVE"/"  19  7*y/.S  VALVE  VE 

249364LL  V-S", 1 0X, " 1 6=TOTAL  HEADLOEE"/”20  *=PIEZ  HD  EELOV  VALVE", ?X, ” 1 7 
2499 04= INEPTI AL  TERM"/"21  9  =  D/S  VALVE  WELL  V-S", I 0X, "13*V- S  D1 FFE°ENT1 A 

2S0004L"/"£2  t’NITS  FOP  AEOVE  PARI AELES  APE*  ”/"23  MI NVTES* 1 ”, 8x, "FT/MI N 
250104= 1  l"/"24  NONE®  2, 3, 4, 10", 4x, "FT/SEC* 1 2"/"25  CFS=5",  I 2x, "FT® I  3,  14, 

250204 15, 16,  17,  13"/"26  FT- DATVM=  6, 7,8,9") 

25030  VRI TEC  2,  1 4840)  CM, 1*1, IS) 

25040  vpl TEC  2,  1 43  60) 

25050  PETf’PN 
25060  END 


52 


